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ABSTRACT 

The Drinciole of pulse modulation used in case of coherent scatter radars 
(MST radars) and incoherent scatter radars (IS radars) is f1r ^ t disc ^ s ® ed - 
rent detection and the corresponding system configuration is delineated. Antenna 
requirements and design are outlined and the phase-coherent transmitter/ receiver 
system is described. Transmit/receive duplexers, transmitters, receivers, q 
drature detectors are explained. The radar controller, integrator, decoder and 
correlator design as well as the data transfer and the control and monitoring by 
the host computer is delineated. Typical operation parameters of some well-known 
radars are finally summarised. 

1. INTRODUCTION TO PRINCIPLES OF THE ATMOSPHERIC RADAR TECHNIQUE 

This tutorial lecture note aims to give a general summary and overview on 
the MST radar and the incoherent scatter radar technique. With the MST ra ars 
the mesosphere, stratosphere and the troposphere are investigated,’ ^ereas w th 
incoherent scatter (IS) radars the ionosphere Is investigated. It Is not in- 
tended to present here a review for experts working In this f i ® ld -. , ° 

give an overview to newcomers to allow a basic introduction for 
joining this research field. For more specific technical descriptions of the 
radar technique In general and the atmospheric radar technique ™ 
the reader is referred to other books and articles, such as those by SKOLNIK 

( 1970 ), HARDY (1972), WILSON and MILLER (1972). BATTAN t^n^of MST 

STRAUCH (1983) DOVIAK and 2RNIC (1984). Very specified descriptions of MST 
rada^technioues caS be found in Handbooks for MAP (BOWHILL and EDWARDS 1983. 
1984 1986), an explicit and useful explanation of the MST and Incoherent scat- 
ter MU radar was published by FUKAO et al. (1985a, b). and a review of UHF and 
VHF radar techniques for atmospheric research and wind profiler applications was 
recently prepared by ROTTGER and LARSEN (1989). Good overviews on the ^coherent 
scatter radar technique are given by EVANS (1969). BARON (1977) and ^ Q FORS 
(1977). Since Important features of the MST and IS radar hardware result from 
certain principles of the basic radar technique, the scattering mechanisms and 
the data acquisition procedures, these will be briefly outlined in the course of 
this lecture note, which, however, mainly aims towards the description of radar 
instrumentation. Further tutorial introductions to the ^ f *" d 

IS radars in atmospheric and Ionospheric research can be found in the other ar 
tides published in this volume of the Handbook for MAP. 

We will describe here the very basic technical principles, which are pro- 
found In all radar applications but are usually adjusted according to specific 
requirements of certain scientific experiments and operational 
Some of those are outlined in the following chapters and for more detal 
qlvenln this tutorial or the in mentioned reviews SaSr 

relevant literature summarized in the reference list at the end of this p pe . 
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1.1. Basics of Pulsed Doppler Radars 

Usually the HST radars as well as the Incoherent scatter radars apply the 
conventional pulse modulation technique, l.e. a short radar pulse Is transmitted 
as shown In the sketch of Fig. 1, and the backscattered radar echo from a range 
r (or altitude z) is received after the time t. Sampling the received echoes as 
function of time then allows to evaluate the echoes from different ranges 
r =ct/2, where c Is the propagation velocity of the radar signal, namely the 
speed of light for the radars operating In the VHF and UHF bands. 



fo 


Figure 1. Principle of a pulsed Doppler radar: 

A transmitted radar pulse is scattered by some refractive 
index irregularity at the ranges r. The backscattered radar 
signal is received after the time of flight t from the 
ranges r = ct/2, where c is the speed of light. Usually the 
power or Doppler spectrum (as shown in the lower Insert) is 
computed for signals received in certain range gates and the 
basic parameters total power P, Doppler shift fo and the 
spectrum width a are deduced. In addition further useful 
parameters can be determined from the particular shapes of 
Doppler spectra. 

The procedure of pulse radar is described in some more detail in Fig. 2. 
Let a pulsed electromagnetic wave be transmitted at the time Ti. The pulse 
duration of this radar signal shall be 5tt . For simplification the pulse shape 
Is assumed to be rectangular, but in real applications it may be a smoothed 
trapezoid or triangle or Gaussian shaped. In a nondispersive propagation medium 
(at the high frequencies used in the MST and IS radar application the refractive 
index is very close to one such that dispersion does not have to be considered) 
the pulse travels with the speed of light c and reaches the range r* after the 
time ti = r«/c. A target at r» can scatter or reflect the radar signal in some 
directions. A small fraction returns to the location of the transmitter, where 
this radar echo will be received after the time ti 1 = 2ti = 2r*/c. This yields 
the basic relation r = ct/2, which allows determination of the range of any 
radar target by measuring the round-trip time t. This relation holds for mono- 
static radars (transmitter and receiver are at almost the same location). For 
bistatic radars (with the receiver separated horizontally from the transmitter 
by a distance comparable to or larger than the ranges to the target) a modified 
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Figure 2. Range-time diagram explaining the main features 
of pulse transmission and scattering from a volume target 
at the range r and the reception after a time t (see text 
for the explanation of the other parameters). 

expression has to be applied. Bistatic or multistatic operation will not be 
discussed in detail here, but is just mentioned for completeness (see Evans, 
1969, for more details). 

Since the transmitted pulse has a finite duration 6tt , Its trailing edge 
will reach the range r* at a time ti+6tt, and reach the receiver at 2ti+6tt = 
ti'+6tt. If the radar target is a thin reflecting surface (called discrete, 
single or hard target scattering), the shape of the received pulse Is the re- 
plica of the transmitted pulse. Now assume that many scatterers fill all ranges 
along the radar beam (called soft target or volume scattering). It is then 
obvious from Fig. 2 that echoes from the range between r*-6r/2 and r* reach the 
receiver simultaneously at ti ’ , and those from r a to r a +6r/2 are received at 
ti’+6tt. The pulse of duration 5tt , thus, at one time illuminates a volume at r a 
extended along a range 5r = c-6tt/2. This is the range gate or range cell from 
which the radar echoes are received. Fig. 2 shows that most of the echo power 
results from the range r a , and minimum power is received from r*±6r/2. Thus, the 
resulting range weighting function of the single range gate centered around r a 
Is a triangle. 

Because of the finite receiver bandwidth (usually matched to the bandwidth 
of the transmitted pulse) the receiver gate has a finite width In time. The 
receiver pulse response (given by about the Inverse of the receiver bandwidth) 
distorts the received pulse by smoothing Its leading and trailing edges and thus 
delaying it by about the receiver response time. The transmitter pulse length 
5tt, to which the receiver response time is matched, determines the range reso- 
lution 5r. The instrumental distortions of the radar pulse are not depicted in 
Fig. 2 to maintain the clearness of survey, but they eventually have to be 
considered in the final analysis of the radar data. 

In radar applications short pulses are normally transmitted periodically, 
i.e. the n-th pulse follows the (n-l)-th pulse after a specified time. For 
convenience this time period is set here to be a multiple (K>1) of 6tt. This 
time (Tn-Tn-O is called the interpulse period Tt PP , IPP, Ti or just T. Its 
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Inverse is called the pulse repetition frequency f p rf = l/Ti PP . The off-on-ratlo 
of the transmitter Ti PP /6tt-l determines approximately the range from which 
radar echoes can be unambiguously received (in units of range resolution). It Is 
more customary, however, to use the ratio d = 6tt/Ti PP , which is called duty 
cycle or duty factor. The average transmitter power Pa, to be averaged over 
(more than) one Interpulse period, is the product of the duty cycle and the 
transmitter pulse peak power P p , i.e. Pa = d-P P . In Table 1 the basic technical 
terms used In radar applications are summarised. 

TABLE 1 

DEFINITION OF TECHNICAL TERMS 
USED IN ATMOSPHERIC RADAR EXPERIMENTS: 



T = interpulse period (IPP), also Ti PP 
1/T = pulse repetition frequency (PRF) 

6tt = transmitter pulse length (duration) 

P p = transmitter peak power 

5tt/T = d = duty cycle 

dP P = average transmitter power 

r = cto/2 = one-way distance to the radar target 


Assume that the radar echo power is due to volume scatter and that isotro- 
pic scatterers totally fill the radar beam. Then the received radar echo power 
P* is given by the radar equation: 

A Pp 5r r) 

Pa = , (1) 

8nr 2 

where A is the effective antenna area and q Is the radar reflectivity. 

If the radar echo is due to reflection from a large surface of a refractive 
index discontinuity, which is stratified perpendicular to the radar wave propa- 
gation, the received radar echo power Pr is given by: 

Pp A 2 ] 9 | 2 

Pr = (2) 

4fio 2 r 2 

where y is the amplitude reflection coefficient of the surface. For MST radar 
applications the reflection coefficient can be very small, resulting in partial 
reflection. In many cases of radar observations scattering takes place, and 
particularly in Ionosphere observations with IS radars the only mechanism is 
(Incoherent) scattering. The partial reflection mechanism is difficult to be 
distinguished from anisotropic scatter, which both have a pronounced aspect 
sensitivity (dependence of echo power on beam direction) observed with the long 
wavelength MST radars. We refer to detailed explanations of these effects, which 
are for instance outlined in other articles of this book (e.g., HOCKING, 1989). 
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As compared to radar echoes from single hard targets (e.g. airplanes), 
where the echo power is proportional to r -4 , P* and Pr are proportional to r~ z . 
This is simply explained by the fact that the volume or the partially reflecting 
surface or layer of the atmosphere, which is illuminated by the radar beam, Is 
not constant but increases with the square of the range. We particularly have to 
note the direct dependency of the received signal power on the transmitted peak 
power Pp and the antenna area A, which means that both these quantities should 
be optimised. For the reason that the sensitivity of the radars Is directly 
determined by these two parameters, transmitter power and antenna aperture, the 
quantity “power- aperture product P-A" is used often as the main Indicator of the 
radar capability or sensitivity and it allows a better comparison of measure- 
ments done at different radars. Also the term slgnal-to-noise ratio (SNR), which 
is the ratio of signal power P»,r and noise power P n , Is used frequently to 
determine the echo signal strength although this can be misleading because the 
noise level varies particularly with frequency and antenna look angle, with 
receiver front-end sensitivity and with loss rates of the antenna and feed 
systems. 

Because in normal radar operations the pulse repetition frequency is kept 
constant, i.e. the transmitted pulse train is periodic, range-aliasing may 
occur. This ambiguity is depicted in Fig. 2, At time t 2 1 an echo of the pulse 
at T 2 is received from range r» , and an echo of the pulse at Ti Is received from 
range rt>. Of course higher order range-aliasing can occur from ranges rn = 
c(t+(n-1 )Tipp )/2 . Because these echoes return from separate scatter volumes, the 
echo signals are uncorrelated but still their power accumulates in the same 
receiver range gate. If no special arrangements (e.g. pulse-coding, frequency 
changes or non-periodic T* pp ) are being made, the maximum unambiguous range is 
rmmx = c-Tipp/2. The minimum range r»in obviously is given by the pulse duration 
5tt plus some instrumentally entailed transition time between transmission and 
reception. 

Assume that a bulk motion carries the scatterers or reflectors In the 
volume at range r. Because of the Doppler effect, the rate of change of phase 0 
of the returned signal is then d0/dt = 4tt/Do dr/dt , where Ho is the radar wave- 
length. When V’ is the (radial) velocity in direction of the radar signal path, 
V’ = dr/dt. The phase change d0/dt is the angular Doppler frequency Qo = 2nfo f 
which yields fo = -2Wrio. Since the radar signal Is pulsed at a frequency fprt, 
i.e. the radar echo Is sampled at a rate T* pp , this yields the maximum Doppler 
frequency to be resolved by pulse- to- pulse analysis (Nyquist frequency): fo»»x = 
fprf/2 = 1/2T 4 pp . This corresponds to a maximum radial velocity V’««x = 
Oofoma x /2 and V’mx = n 0 c/8r«ax . V’»»x, as defined here, is much larger than any 
realistic velocity observed in the lower and middle atmosphere. 

In cases of large Doppler shift and/or large spectrum widths, i.e. the 
Inverse of the signal correlation time, which are observed for incoherent scat- 
ter echoes from the Ionospheric E- and F-reglon, the maximum frequency cannot be 
resolved by the pulse- to- pulse technique and other means have to be applied. The 
difference of "coherent signals" observed in MST radar applications and ’’incohe- 
rent signals" observed in IS radar applications is sketched in Fig. 3. For short 
signal correlation times (shorter than the Interpulse period) all samples for 
the autocorrelation function and the spectrum analysis have to be taken during a 
single Interpulse period (incoherent signal). For long correlation times (longer 
than the interpulse period) these samples need to be taken during several Inter- 
pulse periods (coherent signals), which is called "pulse-to-pulse" technique. 
Details of the corresponding Instrumental configurations are discussed in the 
following Chapter 1.2 and of the data acquisition in Chapter 5. 
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Figure 3. Range-time-ampl Itude diagrams for coherent (underspread) and 
Incoherent (overspread) signals. The dots mark the sample points, which 
are used to compute autocorrelation functions ACF. The tables symbolize 
that for coherent signals coherent integrations (11,12,...) are performed 
for fixed range gates over several interpulse periods, then the ACFs are 
computed which are finally accumulated (XACF). This is called "pulse-to- 
pulse technique ,, . For Incoherent signals the ACFs are computed within 
every interpulse period (ACF1 , ACF2,...) and are accumulated (XACF). 

1.2. An Overview of Particular Radar Systems 

Radars operating in the extended frequency range from MF (medium frequency) 
to UHF (ultra-high frequency) are used to Investigate the structure and dynamics 
of the troposphere, stratosphere, mesosphere and thermosphere/ionosphere. In 
Table 2 different kinds of radars and their basic technical parameters are 
summarized. The MF radars, making use of partial reflection from electron den- 
sity irregularities, are particularly applied to measure electron density pro- 
files of the Ionospheric 0-region and lower E-region as well as the Horizontal 
wind velocity in this altitude range, comprising the mesosphere and lower ther- 
mosphere. The HF radars (conventionally known as ionosondes) and the ionospheric 
irregularity scatter radars (somewhat ambiguously also called "coherent radars") 
are used to study total reflections from the ionosphere and scattering from E- 
and F-reglon plasma irregularities. The meteor radars make use of echoes retur- 
ned from meteor trails to measure wind velocities in the upper mesosphere and 
lower thermosphere. The mesosphere-stratosphere- troposphere (MST/ST) radars, 
detecting echoes from turbulence- Induced inhomogeneities of the radio refractive 
index, are applied to investigate winds, waves, turbulence and stability in the 
Indicated altitude regions. Incoherent scatter (also called Thomson scatter) 
radars, making use of scatter from free electrons, are applied to study the 
Ionosphere and thermosphere. IS radar echoes from the 0- region also have long 
coherence times like MST radar echoes. We concentrate here on the basic prin- 
ciples and applications of the Incoherent scatter (IS) and the MST radars, but 
we note that many principles apply also to the other kinds of radars. 
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TABLE 2 

RADAR METHODS FOR INVESTIGATIONS OF THE LONER AND MIDDLE ATMOSPHERE 
AND THE THERMOSPHERE/ IONOSPHERE 

Typical operation parameters (approximate} 


Radar 

Method 

Frequency 

Range 

Wavelength 
A in « 

Average 
Power 
in IN 

Antenna 
Dimension 
in 1) 

Height 

Region 

MF Radar 

HF-HF 

150-50 

0.01-1 

HO 

H,l 

HF Radar* 

HF 

300-10 

0.01-5 

0.5-1 

Th/Io 

Coherent Radar 4 

HF-VHF 

30-1 

0.1-1 

5-50 

Th/Io 

Meteor Radar 

HF-VHF 

10-6 

0.1-10 

2*10 

M,IT 

MST Radar 

VHF 

6-7 

1-100 

5-50 

M,S,T 

Incoherent 
Scatter Radar 

VHF -UHF 

1.4*0.25 

100-300 

100-300 

H,LT 

ST Radar 

VHF -SHF 

6-0.1 

1-500 

10-500 

S,T 


Hf = 0. 3-3.0 mhz 
HF * 3.0-30 MHz 
VHF = 30-300 mi 
UHF = 300-3000 MHz 
SHF = 3-30 GHz 
* = Ionosonde 


H = Mesosphere 

S s Stratosphere 

T = Troposphere 

LT = Lower Thermosphere 

Th/Io = Thenosphere/Ionosphere 

♦ = Irregularity Scatter 


The MST radars operate In the lower VHF band around 50 MHz, corresponding 
to wavelengths around 6 m. Since quasi-vertical antenna beam directions are 
used, ranges are roughly equal to altitudes. For MST radar observations of the 
lower and middle atmosphere the range limits r««x are between 10-20 km and 100 
km. This yields typical pulse repetition frequencies between 10 kHz and 1 kHz. 
Altitude resolutions from about 1 km down to at least 100 m are required to 
resolve typical vertical scales in the troposphere, stratosphere and the meso- 
sphere, This corresponds to pulse lengths of about 1-10 us. Thus, typical duty 
cycles are between about 10- 1 and 10 -3 . Longer coded pulses are frequently used, 
which is discussed in Chapter 4.5. 

MST radars make use of scattering and reflection from variations of humi- 
dity, temperature and electron density, Induced by turbulence In the lower and 
middle atmosphere. Essentially, MST radars can observe: the 3-dimensiona I wind 
vector, atmospheric reflectivity and stability, and morphology of turbulence and 
waves. The continuous measurements with MST radars offer very good quality and 
quantity middle atmosphere observations of wind velocities. MST radars operate 
at frequencies around 50 MHz, and therefore are also called VHF radars (VHF = 
very high frequency band between 30 MHz and 300 MHz). Higher frequency radar* 
mostly cover only the troposphere and stratosphere. Typical peak powers of VHF 
radars are between 1 kW and 1 MW. Range resolutions down to about 100 m and tla* 
resolutions down to some ten second are possible. The antenna arrays with typl- 
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resolutions down to some ten second are possible. The antenna arrays with typi- 
cal dimensions of 1000 m 2 to some 10000 m 2 point close to the zenith direction. 
Coherent detection, digital control and data acquisition are mandatory as is 
described In Chapter 4. 

There are basically two methods which are applied by MST radars. In parti- 
cular to measure velocities, as is sketched In Fig. 4. One method uses a narrow 
radar beam pointed into various directions and measures the Doppler shift of 
echoes scattered from irregularities. This method Is usually called the "Doppler 
method" and for this reason these radars are also called "Doppler radars". The 
other method uses three or more spaced antennas and the received echoes are 
cross-correlated to determine the drift speed of irregularities and is called 
"spaced antenna or SA method". Since the irregularities are usually moving with 
the wind velocity, both methods are capable to measure the wind velocity. Al- 
though both methods are basing on the same physical mechanism, in praxis the 
technical implementation may favour one or the other of these methods. The 
spaced antenna method can also be applied in the spatial domain radar interfero- 
meter mode, which is advantageous for studying the structure of the scattering/- 
reflecting irregularities. In all these applications a phase-coherent radar 
system is applied, it is even required for the Doppler and the interferometer 
method . 


3-DIM VELOCITY MEASUREMENTS WITH VHF-RADAR 



Figure 4. The two principle methods of three-dimensional velocity measurements 
with monostatic mesosphere-stratosphere-troposphere (MST) radars. These 
methods are the Doppler method using oblique beams to deduce the horizontal 
and vertical velocity from the measured radial velocity, and the spaced 
antenna method apply-ing drifting pattern measurements to deduce the horizon- 
tal wind velocity. 
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The smaller versions of this Doppler radar category, the stratosphere- 
troposphere (ST) radars operate according to the same principle like the MST 
radars, but frequently operate on higher frequencies in the UHF band with smal- 
ler powers and antenna apertures of several ten square meters only. They are not 
capable to detect echoes from the mesosphere since either their power- aperture 
product Is too small, yielding too low a sensitivity, or their frequency is so 
high that their wavelength is In the Inertial subrange of refractive index 
variations such that the scatter echoes are not detectable. These radars apply 
the Doppler method to measure velocities and are mostly very capable to investi- 
gate the lower atmosphere (particularly the troposphere and the boundary layer) 
with a very good height resolution of some 10 meters or better. 

The state of the art of the ST radars and their suitable applicability to 
measure continuously the wind profiles in the troposphere and the lower strato- 
sphere has attracted wide attention in the meteorological community. As a conse- 
quence such ST radars are included Into meteorological research networks and 
particular Instrumentation Is being designed to allow long-term routine opera- 
tions of these so-called wind profilers. Since this Is a development by in- 
dustrial companies, we will not discuss these particular designs In the context 
of this tutorial. A short summary of system specifications of wind profilers, 
however, can be found in Tables 6 and 7 in Chapter 6. 



Figure 5. The principle of tri-static incoherent scatter 

(IS) radar measurements of the auroral ionosphere by the EISCAT 
UHF system with transmitter and receiver in Tromsd, Norway, 
and remote receivers in Kiruna, Sweden, and Sodankyla, Finland, 

The incoherent scatter radars operate at VHF and UHF frequencies with a 
much larger power- aperture product than the MST or ST radars. Mostly these 
radars have large dish antennas or phased arrays occasionally with diameters of 
several ten to hundred meters. Their peak powers are usually larger than 1 MW, 
their duty cycle can be up to 12 percent and their fastest interpulse period is 
mostly not smaller than 1 ms. Pulses as long as 1 ms are used, but frequently 
amplitude- or phase-coded. Typical altitude resolutions are 0.3 -50 km and time 
resolutions of some minute to several 10 minutes are achieved. These differences 
result from the fact that because of the very small volume reflectivity of In- 
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coherent scatter, the signal-to-noise ratios are mostly well below unity, where- 
as the SNR Is usually much larger than 1 in the MST radar case. The incoherent 
scatter radars particularly measure electron density, electron temperature, ion 
temperature and ion velocity in the ionosphere. Fig. 5 shows the tri-static 
Incoherent scatter UHF-radar system of the EISCAT Scientific Association which 
Is operated In northern Scandinavia to study the auroral ionosphere. The tri- 
static radar configuration is most useful to measure the three-dimensional 
velocity with the best achievable accuracy. All the IS radars were also used for 
studies of the stratosphere and the troposphere, although some Instrumental 
constraints had to be considered such as ground clutter, receiver recovery and 
dynamic range limitations. 

A list of these radars is found in Table 9 at the end of this article. 

1.3. The Principle of Phase-Coherent Radar Systems 

Incoherent or Thomson scatter is caused by fluctuations in the radio re- 
fractive index at the radar Bragg wavelength resulting from thermal motions of 
free electrons in the presence of Ions In the ionosphere. Due to these random 
motions the Incoherent scatter signal is widened in frequency and it needs to be 
sampled fast enough to resolve the complete Doppler spectrum. In the lower 
ionosphere collisions between ions and neutrals Impose the fluctuations of the 
neutral atmosphere on the plasma and consequently on the electrons. Since these 
fluctuations, although they can be strong in amplitude, do not contain such high 
frequency components as the thermal motions, the incoherent scatter signal 
displays a narrower spectrum, l.e. a longer coherence time in the lower iono- 
sphere. Eventually, the structure of the plasma In the mesosphere (D-reglon) 
will be completely governed by the fluctuations due to neutral air turbulence. 
Instead of Incoherent scatter this latter process is called coherent or turbu- 
lence scatter. Usually the signal spectrum due to turbulence scatter is much 
narrower than that of Incoherent scatter. In addition to the widening of these 
signal spectra, they are shifted in frequency If there is a bulk radial velocity 
of the scattering medium. We would like to explain here that both these scatter- 
ing processes are being studied with phase-coherent radar systems. 

Knowing that fluctuating velocities and radial velocities with quasi-verti- 
cal MST radar beams do not exceed several 10 ms -1 , the Doppler frequency will 
barely exceed 10 Hz for MST radars operating at VHF and will be ;just an order of 
magnitude larger for UHF. Applying fprf > 1 kHz, the MST/ST radar echo will be 
oversampled, l.e. its phase and amplitude does vary little from pulse to pulse 
as it was sketched in Fig. 3. This Is called a coherent radar echo, in contrast 
to an Incoherent radar echo, which, because of the much shorter correlation or 
coherence time of the ionospheric scatter medium, randomly changes phase and 
amplitude from one pulse to the next. Also the velocities in the ionosphere 
measured with incoherent scatter radars is usually much larger than the veloci- 
ties in the middle and lower atmosphere, which additionally causes a more rapid 
change of the phase of the incoherent scatter signals. One does make efficient 
use of the characteristic coherency of MST radar echoes to improve the data 
acquisition procedures, which will be discussed in Chapter 4. These differences 
of the signal coherence times also determine the principle way of the data 
acquisition structure and instrumentation. 

The differences between "coherent signals" detected by the MST/ST radars, 
used to study the lower and middle atmosphere, and "incoherent signals" detected 
by the Incoherent scatter radars, used to study the ionosphere, is elucidated by 
the schematics of Figs. 6a and 6b. Since the coherence time, i.e. the inverse of 
the width of the Doppler spectrum, generally Increases from the lower to the 
middle and to the upper atmosphere (ionosphere) due the increasing variations of 
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MST-RADAR 



Figure 6a. Block diagram showing the basic configuration of an 
MST radar system (data acquisition for coherent signals). 


IS-RADAR 



Figure 6b. Block diagram showing the basic configuration of an 
incoherent scatter radar system (data acquisition for Incohe- 
rent signals). 
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the scattering medium. Whereas In the lower atmosphere the coherence time (de- 
termined principally by turbulent fluctuations) is generally much longer than a 
typical Interpulse period of an MST/VHF radar of say 1 ms. It Is generally much 
shorter (given by thermal fluctuations) than an Interpulse period of Incoherent 
scatter radars. This means that coherent Integration and the autocorrelation 
functions (ACFs) or Doppler spectra are computed on the pulse- to- pulse basis for 
coherent signals ("pulse-to-pulse technique”). On the other hand, the ACFs have 
to be Immediately computed without coherent Integration for each single Inter- 
pulse period using single long or coded pulses In the case of incoherent sig- 
nals. The ACFs (from which the spectra are computed) are then postintegrated. 
The definition of a coherent signal, which we use here. Is consequently that the 
coherence time Is much longer than the Interpulse period. The Incoherent signal 
Is defined by a coherence time which Is much shorter than the Interpulse period. 
At some places In the early Ionosphere literature the terms under-spread and 
over-spread signals are used Instead of coherent and Incoherent signals, respec- 
tively. Both kinds of signals anyhow result from the scattering process being 
coherent In the former and incoherent In the latter case, and both are studied 
with phase-coherent radar Instrumentation. 

We notice that the Instrumental technique as well as the data acquisition 
and analysis of the MST radars as well as the IS radars are basically similar, 
they are well developed, elaborated and fairly mature, although further refine- 
ment Is always going on. We now will briefly outline the two only striking 
differences of the Instrumental design of these two main radar categories, which 
result from the coherent and incoherent signal properties and later will discuss 
particular parts of the Instrumental hardware and their principles. 

In order to allow the measurement of all signal parameters, particularly 
the Doppler spectrum, all MST radars and IS radars are phase-coherent. This 
means that the receiver (RX) and the transmitter (TX) are controlled by or 
phase— locked to the same main oscillator (see Figs. 6a and 6b). The transmitter 
radiates through the antenna the electromagnetic signal S, which propagates to 
the radar volume In the atmosphere or the Ionosphere. From there a small frac- 
tion of the electromagnetic energy Is backscattered to the antenna (which can be 
a separate or the same antenna, provided that a transmit/receive duplexer can be 
used). External noise, received by the same antenna and Internal receiver noise 
adds to the atmospheric/ionospheric echo and the total signal plus noise C’ Is 
amplified in the receiver and mixed to base band by the same oscillator signal 
S°, which Is used to control the transmitter. The resulting base band signal C 
Is' converted to a digital series in the analog-to-dlgltal converter (ADC), which 
is controlled by a series of sample pulses from the radar controller. The radar 
controller also generates the transmitter modulation and other control pulses. 

Following the ADC, the data acquisition procedures are different for MST- 
and IS-radars as delineated In the Figs. 6a and 6b. This difference results from 
the different coherence times as described earlier. Because of the long cohe- 
rence times of echoes from the mesosphere, stratosphere and troposphere, the 
digital data are usually coherently added In the Integrator of an MST radar 
(Fig. 6a), followed by a decoding procedure If the transmitter pulses are coded 
(see Chapter 4 for details). These raw. but pre-lntegrated and/or decoded, data 
are either directly dumped Into the host computer or analyzed In terms of auto- 
correlation functions or Doppler (Fourier) spectra. The host computer Is used to 
store the raw data together with other operating system parameters and the time, 
on file or tape, perform an on-line qulcklook analysis or some further prepro- 
cessing such as computation or post-integration of ACFs or spectra. The host 
computer also supervises the entire system by loading and starting the radar 
controller etc. as well as checking the system performance and Issuing Inter- 
rupts In case of system malfunctions. In the case of an Incoherent scatter radar 
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(Fig. 6b) first a decoding has to be performed, then usually ACFs are computed 
on-line which are post- integrated and then dumped into the host computer. In 
particular pulse coding schemes also an additional decoding procedure has to be 
applied after the ACF computation. Otherwise the instrumental design of an IS 
radar is principally similar to an MST radar. 

2. RADAR ANTENNAS 

The radiation pattern of MST- and IS-radar antennas has to be carefully 
designed in order to optimize the wanted atmospheric/ionospheric scatter signal 
as compared to interfering external influences. These are summarized in the 
sketch of Fig. 7. It is evident that all the unwanted components, - ground 
clutter, sea clutter, airplane or satellite echoes, radio interference and 
scatter received through antenna sidelobes - , have to be properly eliminated or 
minimized. To achieve a high sensitivity and a reasonable angular resolution the 
antenna gain should be large and thus the antenna beam width small. The cosmic 
noise level picked up by the antenna is unavoidable in the low VHF band. In the 
high VHF and in the UHF band the noise level is not given by the sky noise but 
by the system noise. This can be minimized by optimum design of the receiver 
front-end amplifier and adapted data processing. 

The antennas of MST- and IS-radars usually consist of an ensemble of single 
elements which are phase-coherent ly combined to a large antenna array or are 
large dish antennas. Usually the dimensions of MST radar antennas are more than 
several ten radar wavelengths, their gain is mostly exceeding 20 dB, correspond- 
ing to a beam width of less than 10 degree. The IS radar antennas usually have 
antenna gains larger than 40 dB In order to detect the weak incoherent scatter 
signals. The beam directions can be mechanically or electrically positioned 


Cosmic Noise 



Figure 7. Schematic view of an atmospheric radar in an interference and clutter 
environment. The radar antenna is depicted here to consist of a phased array 
of Yagis, but could principally be any other kind of antenna system. 
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Figure 8. Near-zenith pointing 
antenna beams of a 403 MHz 
phased array antenna of an 
ST radar (after FRISCH 
et al., 1986). 


(e.g., Fig. 8) in order to observe spatial variations of the scatter process and 
to measure three-dimensional velocity components. The principle of measurements 
of the velocity components is explained by the schematics of Fig. 9. At a zenith 
angle 6 in the direction of the wave vector k the radial velocity v = u’+ w* is 
measured, which consists of the projections w* of the vertical velocity w and 
the projection u’ of the horizontal velocity u. The combination of the veloci- 
ties measured In the directions 5 and -5 yields the horizontal and the vertical 
velocity component. Depending on the requirements to measure all velocity compo- 
nents and their divergence as well as vorticity, more beam directions are ne- 
cessary. Multi-static radars (e.g., Fig. 5) can also be used to measure the 
velocity components. The discussion of these possibilities and requirements is 
out of the scope of this lecture note but can be found elsewhere (see ROTTGER 
and LARSEN, 1989, who discuss the applications of MST radars and give further 
references) . 



V1-V2 
2 sin 8 


w = 


V 1 4 v 2 
2 cos 8 


Figure 9. Principle geometry of two antenna beam directions koi and ko 2 at 
zenith angles 6 and -6 for the deduction of the horizontal velocity compo- 
nent u and the vertical velocity component w. The measured radial velocity 
components vi and V 2 consist of the projections u* and w’ of the u and w wind 
velocity components, respectively. 
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We now will briefly describe the basic antenna theory in order to allow an 
understanding of the antenna functions and the design optimization and then will 
describe some typical antenna systems used in MST radar applications. 

2.1. Some Basic Antenna Parameters 

The main parameter determining the antenna gain and the corresponding beam 
width is the size or the area of the antenna, which is also called the antenna 
aperture. The antenna aperture can be either regarded as a continuous radiator 
If it is illuminated by an independent feed system, or it can consist of several 
sub-elements which are electrically fed from a common transmitter source. Both 
approaches can be basically treated by a similar formalism, which we will ex- 
plain for the simple case of a line of antenna dipole elements. Assume that N 
individual antenna elements with equal spacing d are horizontally lined up to 
form a multi-element array (see Fig. 10). In the array far-field (r > (Nd) 2 /n 0 ) 
the antenna polar diagram is 

E(6) = I En (5)exp( 1 ( ------- s1n6 + ♦n)), (3) 

flo 

n=1 

where En(6) Is the pattern of an individual element and ♦n Is a relative phase 
placed on this element. Thus, the polar diagram is just the Fourier transform of 
the spatial aperture distribution. If all the elements have similar En(5) and 
4>n, the polar diagram Is a function sin Nx/sin x with x = ird -sinS/flo . If d < 
flo, there Is only one main lobe at 5 = 0° (for 0 = 0). If d > flo, grating lobes 
at 5g = arcsinOWd) occur, for En(6) = const with amplitudes similar to the 
main lobe. The width of the main lobe is 6b = arcsin (0 o /Nd). In radar applica- 
tions also the two-way beam width 6b/V2 is used. If 6b Is small, it is directly 
proportional to the ratio of the wavelength flo to the aperture dimension Nd. The 
radiation pattern has minima (nulls in theory) at 6n = ± p6b, where n = 1,2,..., 
N/2, and sidelobes occur at 6» = 6n ±6b/2. If equal weighting Wn is applied to 
each of the single elements (e.g., En’(6) = WnEn(5), with Wn = const for all 
n), the first sidelobe closest to the main lobe is suppressed by 13.2 dB plus 
the fall-off of the Individual element pattern En(6) at 6»i . 



Figure 10, Schematic drawing of wave vectors of a plain wave 
radiated under a zenith angle 6 from N isotropic antenna 
elements with spacing d. 
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The antenna pattern, particularly the desired variation of the direction 60 
of the main lobe, can be changed by applying a linearly progressing phase <t>n 
from element to element, which has to be $n = 2itd(n-1 )sin 6 o/n 0 . This beam 
steering should be within reasonable limits of the individual element pattern 
En(6) to avoid undesirable degrading of the antenna radiation. Furthermore the 
effective area of an antenna, which is fixed on the ground, is reduced with 
larger zenith angles because the illumination changes by sin 6 . 

To obtain improved sidelobe suppressions, a tapering of the antenna array 
can be applied by either changing the weighting function Wn (electrical weight- 
ing, i.e. feeding the outer elements of an array with less power than the Inner 
elements), or by using unequal element spacings dn (spatial weighting, i.e. 
applying larger spacings for the outer elements of an antenna array). The price 
one has to pay for the improvement of sidelobe suppression by tapering is always 
a broadening of the main lobe and lowering of the antenna gain G. Using a tri- 
angular weighting, for instance, improves the one-way suppression of the first 
sidelobe to -26 dB, but widens the main lobe by a factor of 1.44 and reduces the 
gain by 25%, as compared to uniform weighting. The respective values for a cos 2 - 
weighting, which is a good approximation to a Gaussian weighting, are -32 dB, 
1.64 and 33%. 

These one-dimensional considerations can be extended easily to a realistic 
two-dimensional radar antenna array by using. Instead of d, the projections d’ 
of all element positions onto n axis elongated in the azimuth direction a. The 
total antenna pattern then can be calculated by (3) for any a and 5. For a real 
radar antenna we also have to consider that the radiation can be into only one 
half-sphere. For an array system this means that reflections from the ground, a 
screen or reflector elements have to be included in the calculations. The an- 
tenna gain G is defined as the ratio of the maximum radiation intensity (in the 
main beam) to the average radiation Intensity (averaged over all 6 and a). For 
an antenna array with reasonable sidelobe suppression it is proportional to the 
antenna area and is given approximately by the inverse of the two orthogonal 
beam widths 6b’ and 6b’’ (in radians): 

4 tt A An 

G = * . 14) 

no 2 6 b '6b” 

The effective area A or aperture of the antenna is the product of the physical 
area of the antenna and the efficiency of its illumination, which for instance 
may be reduced by tapering. It is noted that A does not include the losses of 
the antenna elements and its feed system. 

The considerations of antenna arrays, consisting of several discrete ele- 
ments, can generally be extended also to antennas with continuous aperture 
illumination, such as dish antennas. The aperture A of a phased array or a dish 
antenna is used to calculate the power- aperture product PA, which defines the 
sensitivity of an MST radar. Note that P is the effective average power, which 
is radiated by the antenna, I.e. it Is smaller than the output power of the 
transmitter due to losses in the antenna and feed lines. 

This very brief outline Is useful for estimating some basic values needed 
to plan and design a radar antenna system. The minimum requirement for investi- 
gations of the troposphere and stratosphere is a power aperture product P-A = 
10 e Wm 2 , where P = Pa is the average transmitter power. We assume that the radar 
sensitivity is a function of A k , where 1 < k < 2 depending on the reflected con- 
tribution (see equation (2)). It means that in practice one should favour an 
extension of the antenna area against an increase of the transmitter power. For 
a commonly achievable transmitter power P = 10 3 W, PA = 10® Wm 2 yields the 
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dimension of a circular antenna array to be about Nd % 36 m, and A * 1000 m 2 . 
For Ho = 6 m we can obtain a half-power beamwidth 6b * 9°, and a gain G % 27 dB. 
This estimate assumes that the array is optimally filled with elements. A single 
dipole element over a proper reflector screen has an effective area of about 15 
m 2 (at flo = 6 m). Thus, about 64 dipoles are needed to fill the array. For a 
square array with 8x8 elements the resulting spacing Is d * 0.75 Ho, and a 
grating lobe will not occur. The first nulls are at 5n = ±9°, and the first 
sldelobes at 6* = ±13.5°. To use such an array to measure reliable velocities, 
the antenna beam has at least to be steered to a zenith angle 6o = 9° to place a 
null into zenith direction (minimizes influences of aspect sensitivity). Even 
then a sidelobe at -4.5° is a problem. 


When designing a radar antenna for atmospheric research one has to trade 
between the choices to optimize the effective aperture or to optimize the side- 
lobe suppressions. An optimization of the aperture increases the sensitivity, 
but does not minimize the side lobes. Suppression of sldelobes by tapering 
attenuates undesirable signals which spoil the estimates of reflectivity and 
velocity, but it reduces the antenna gain. Principally, any sidelobe effects or 
their minimization, however, are equivalent to a broadening of the antenna beam. 
This Is generally a minor shortcoming as compared to the tedious procedures to 
remove sidelobe effects from the signal during the data analysis procedure. As 
an example of a well optimised HST radar antenna pattern we show In Fig. 11 the 
antenna diagram In two vertical planes of the mobile SOUSY VHF Radar (from 
CZECHOWSKY et al . , 1984). Note the suppression of near-zenith sldelobes by 25 dB 
as compared to the main lobe and the grating lobe at 40°. 


The radar echo signal Is given by the product of the antenna pattern with 
the spatially varying reflectivity structure of the atmosphere. Thus, Knowing 
the antenna pattern, it should be In principle possible to find the wanted 
signal parameters which however is generally an ambiguous and time-consuming 
inversion procedure. For vertically pointing main beams the sidelobe effects are 
efficiently suppressed if there is anisotropic scattering with maximum aspect 
sensitivity in zenith direction (e.g., for the 50 MHz MST radars). It follows 
that sldelobes in such cases are a minor problem for investigations with verti- 
cal beams. However, they can be crucial for methods applying off-zenith beams 
with 50 MHz MST radars. If a sidelobe is pointing towards the zenith a larger 
power may be received from the vertical than from off-vertical directions when 
there is an aspect sensitivity due to horizontally stratified scatterers or 
partially reflecting layers. 



Figure 11. Computed radia- 
tion pattern of the mobile 
SOUSY VHF radar (53.5 MHz) 
phased array, consisting of 
576 Yagi antennas. This 
antenna system is located on 
Andoya, Norway. (From 
Czechowsky et al., 1984). 


o 
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Also sldelobes at low elevation angles have to be considered since these 
can cause strong echoes from the non-atmospheric targets In the surroundings of 
the radar antenna (particularly the ground clutter seen In Fig. 7). The super- 
position of ground reflected (hatched line in Fig. 7) and direct radiation does 
very effectively reduce the radiation at grazing angles (5 > 85-90°) since the 
antenna elements of a phased array are horizontally polarized and the ground-re- 
flected wave then suffers a phase reversal during reflection. This even can 
suppress low sldelobes of the array pattern which may be regarded as crucial If 
one would not take Into account ground reflections. The location of an MST radar 
antenna array at a flat ground (extending out to several 100 m) may be suffici- 
ent, but a shallow valley generally should be preferred to further suppressing 
the low angle radiation effects. However, high extending targets, such as .radio 
towers or mountains In the close vicinity, will still cause considerable clutter 
echoes, even when optimising the antenna array for low angle radiation suppres- 
sion. For IS radar antennas the criteria of low angle radiation suppression are 
not so stringent, because ground clutter barely comprises a problem since the 
Ionospheric target regions are usually at much larger ranges. However, for IS 
radars In auroral regions care has to be taken to suppress coherent echoes from 
Ionospheric irregularities. 

2.2. Antenna Types and Feed Systems 

MST radars generally operate with quasi-vertical beams, l.e. the zenith 
angles are smaller than about 20°-30°. As shown In Fig. 8 usually zenith angles 
of 15° are chosen. IS radars mostly apply a variety f beam directions, also 
those very close to the horizon. For MST radar operation linear polarization Is 
sufficient, but circular polarization Is needed for Incoherent scatter radars 
since the polarization can change due dispersion In the Ionosphere. Essentially 
four different types of antenna systems are In use: dish antennas, dipole 
arrays, £oax1al-£ol 1 Inear arrays (COCO) and Yagl arrays. In Fig. 12 we show as 
an example the 32-m dish antenna of the EISCAT Incoherent scatter radar 933-MHz 
receiving station In Sweden, which can be pointed in all directions by mechani- 
cal steering. In Fig. 13 the 200m x 200m MST radar antenna of the Poker Flat 
Radar in Alaska Is shown which has only a few fixed pointing directions to the 
vertical and at 15° off-vertical. 



Figure 12. The 32-m receiving 
antenna dish of the Incohe- 
rent scatter radar of the 
EISCAT UHF (933 MHz) radar 
system in Klruna, Sweden. 
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Figure 13. The 200m x 200m coaxial-col 1 inear antenna array of 
the Poker Flat MST radar in Alaska, operating on 50 MHz with 
two off-vertical beam directions (from BALSLEV et al., 1980). 

Dish or cylinder antennas are rarely used for MST radar applications be- 
cause of their large dimensions, but they are quite frequently applied in IS 
radar facilities. Only one or a few elements are applied as primary feed anten- 
nas of the dish antennas. The beam steering is done either by moving the posi- 
tion of the feed antenna (e.g., at the Arecibo Observatory (WOODMAN, 1980)) or 
the entire dish (e.g., at the other UHF IS radars, see Table 9 at the end of 
this article). This has the advantage that no complicated power distributions 
and phasing network has to be used to feed the antenna. The EISCAT VHF radar 
uses a long line feed to illuminate a 120m x 40m dish and allows mechanical 
steering in one plane and side-steering by phasing in the other plane (HAGFORS 
et al., 1982). Because of the limited size of the primary feed antennas of dish 
antennas, the low-angle sidelobe suppression is usually not sufficient, which 
results in strong clutter echoes, particularly when such antennas are also used 
for MST radar applications. 

Usually the antenna systems of MST radars are phased array antennas con- 
sisting of many single elements as can be seen for instance in the lay-out of 
the antenna system of the MU radar in Japan which consists of 475 crossed Yagi 
antennas (Fig. 14). In the case of the MU radar modules of four antenna elements 
are fed by single phase-coherent transmitters, which allows very fast electronic 
beam steering. In many other cases the individual elements or modules are fed 
from a single transmitter through a cascading network of cables, hybrids and 
phase shifters. A relatively simple antenna array is the COCO antenna (coaxial- 
col linear) , which is built in form of an array (see Fig. 15) by using coaxial 
cable as radiating elements (e.g., BALSLEV et al . , 1980). It has the advantage 
that the feeding of elements in one line is lust done by interchanging the inner 
and outer conductor of a coaxial cable every half wavelength. The earliest 
application of the COCO antenna can be found at the Jicamarca VHF radar in Peru 
(see Fig. 16). The outer conductors of the aligned coaxial tube or coaxial cable 
act as col 1 inear dipoles. The feeding is done from the center of a line, which 
may typically consist of 16-48 dipoles. Positioning several of these strings or 
rows in parallel at spacing d < Ho , and feeding these by a suitable matching 
network results in a COCO array (see Fig. 15). The radiation and the loss in a 
coaxial string comprise some natural tapering, having the intrinsic advantage of 
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Figure 14. The lay-out of the 475 
crossed Yagi-antenna elements of 
the 46.5-MHz MU radar in Japan 
(after FUKAO et aU, 1985a, b). 



Figure 15. The coaxial-col 1 inear 
antenna lay-out of the 50-MHz 
Sunset ST radar in Colorado, USA 
(from Green et al., 1975). 



Figure 16. Part of the phased-array 
antenna of the 50-MHz Jicamarca 
radar in Peru and a close view of 
the interchange of Inner and outer 
connector of the coaxial-col 1 inear 
dipole tubes (from OCHS, 1965). 


Figure 17. Yagi antennas and power 
dividers of the SOUSY-VHF-Radar 
in West Germany. 
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suppressing sidelobes in the plane of the string but the disadvantage of degra- 
ding the antennna efficiency. Because of the phase relation along a string Is 
fixed, a beam steering parallel to the COCO string is not possible. Beam steer- 
ing perpendicular to the string rows is achieved by inserting appropriate phase 
delays In the cables feeding the parallel rows. The COCO dipoles have to be 
lined up a quarter wavelength over reflector wires or screens. This antenna type 
is relatively inexpensive, because coaxial cable Is used for antenna elements 
and the matching network Is simple. The successive phasing from one collinear 
element to the next, however, degrades the bandwidth of this antenna type. 
Instead of coaxial cable as radiator, halve-wave dipoles can also be used which 
are fed In a properly adjusted phase to form a collinear antenna. The applica- 
tion of collinear dipole lines limits the steerability of an array, and for this 
reason frequently three antenna arrays are used with three different fixed beam 
directions (e.g., GREEN et al., 1975; BALSLEY et al., 1980). 

Single dipole or Yagl antenna elements (see Fig. 17) are often set up to 
form a phased array (see Fig. 14 and Fig. 18). They are fed by a cascading 
network of open wire or coaxial cable systems, when a single transmitter is ap- 
plied. The cascading Is most appropriately done In 2 n branchings (n = 1.2,...), 
as for Instance shown in Fig. 19. By these means one can also feed parts of the 
antenna array with 1/2, 1/4,... power to provide tapering without dissipating 
power (see for instance Fig. 20). The branching is best done in couplers, power 
dividers or 3dB-hybrids, which prevent RF-power, reflected from a mismatch to 
return to the transmitter or other antenna elements rather than being dissipated 
in the resistor port of the hybrid. These hybrid installations also minimize 
effects due to mutual coupling between the single antenna elements. The coupling 
can be critical if phase control is applied to steer the antenna beam. The phase 
shifting, as shown in Fig. 21, Is usually Inserted close to the final elements, 
which may also be connected to form modules of 4 or more elements. For continu- 
ous beam steering phase shifting has to be done continuously, which mostly is a 
sensible and difficult procedure particularly when narrow antenna beams are 
used. The phasing is more easily done by inserting discrete phase delays in 
steps of 2n/16 with a binary phase shifter, which is shown in Fig. 22. The 
switching can be achieved by only four relays, allowing phase delays of 22.5°, 
45°, 90° and 180° to be Inserted in all possible combinations, which are suffi- 
cient for commonly applied beam width of several degrees. Another way of phase 
shifting can be done with hybrids as is shown in the schematics of the later 
Fig. 32. 
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Figure 18. Lay-out of one 64-Yagi-antenna module of the 52-MHz 
Chung-Li VHF radar In Taiwan (after BROSNAHAN et al., 1984). 
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Figure 19. Power divider and beam-steering phase delay (0°, 90°, 
180° and 270°) of the Chung-Li VHF radar, allowing five beam 
positions: vertical (6 = 0°) and north, east, south and west 
at 5 = 16.7° zenith angle (from BROSNAHAN et al., 1984). 



Figure 20. One quarter of the power divider system of the mobile 
SOUSY VHF radar antenna, allowing for tapering of tne antenna 
elements by staggered power dividing hybrids (from CZECHOWSKV 
et al . , 1984) . 















76 


Antenna Feed Distribution System 



Figure 21. Schematics of a typical antenna feed distribution system 
using hybrids for power division. The two output ports of a hybrid 
are phase shifted by 90° with respect to each other. This has to 
be compensated by the phase shifters, which are also used to insert 
certain phase shifts for beam steering. The power, which Is reflec- 
ted from the antennas due to an impedance mismatch, is attenuated 
in the resistor port of the hybrids and does not return through the 
transmitter-receiver port. 


Binary Phase Shifter 



Figure 22. Principle of a binary phase shifter, allowing all the combi- 
binatlons of phases: 

Phase shift: tt/8 2n/8 3n/8 4n/8 5n/8 6n/8 Tn/8 n 

Phase delay A: 10101010 

(1=on) B: 0 1 1 0 0 1 1 0 

(0=off ) C: 0 0 0 1 1 1 1 0 

D: 00000001 

The phases between n and 15/8n are achieved by keeping the phase 
delay 0 in the on-position. 
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The advantage of Yagl against dipole antenna elements Is that no ground 
screen Is needed because of the Yagl reflector. The multi-element structure of a 
single Yagl allows for a higher gain (Improving the filling factor of an array) 
and a negligible coupling (<-25 dB) between adjacent Vagi antennas In an array. 
Mostly Yagi antennas can be constructed In such a way that the bandwidth is 
several Megahertz. The bandwidth limiting factors In a Yagl system essentially 
are the phase shifters. The losses of a Yagl system are also considerably lower 
than those of a COCO antenna. However, such a Yagl system is obviously more 
expensive than a COCO system. The feeding of an array system can be from one 
transmitter, but also sub-modules can be fed separately by several phase-con- 
trolled transmitters (FUKAO et al., 1980). The transmitter phase control can 
even be used to steer the antenna beam. This, however, needs a similar phase 
control of the receiver channels. 

3. TRANSMIT-RECEIVE SYSTEM 

Separate antenna arrays can be used for the transmission and for the recep- 
tion mode as shown in the graph of Fig. 23. Such an antenna combination was 
first used with the SOUSY-VHF-Radar for MST measurements with the beam-swinging 
as well as the spaced-antenna mode, which was added at a later time. For a newly 
designed radar a combination of the two measurement modes using one antenna for 
transmission and reception is useful as shown in Fig. 24. While separate antenna 
arrays allow for sufficient decoupling of the receiver from the transmitter, a 
fast and highly Insulating transmit-recelve duplexer (ATR = antenna-transmit- 
recelve switch or duplexer in Fig. 24) has to quickly switch the antenna from 
the transmitter to the receiver and vice versa, if only one antenna Is used. The 
principle of a duplexer is outlined in Fig. 25. During the transmission phase 
the pin-diodes are short-circuited and the power at the hybrid output ports is 
reflected back into the hybrid from where it exits through the antenna port to 
the antenna. The small leakage through the diodes Is combined in phase in the 
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Antenna I 196Yagis 
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(Doppler velocity) 
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Antenna H , 3*32 Yogis 
fixed vertical beams 
(drift velocity) 

gain *3* 22 dB area *3 * 395 m* 


Figure 23. Antenna systems of the SOUSY-VHF-Radar In W. Germany consisting of 
phased arrays of Yagi-antenna elements. The complete system is designed for 
beam steering in the Doppler mode with antenna I and in the spaced antenna 
mode with antenna I and antenna II. For the latter purpose the antenna I Is 
pointed vertical and the separated three subunits of 32 Yagls of antenna II, 
which are forming a vertical beam also, are used for reception. 
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Figure 24. Radar system using three separate receiver-transmitter-antenna 
modules, which can be operated phase-synchronized to allow beam steering as 
well as spaced antenna measurements with the same full antenna array. The 
system is laid out for range multiplexing (RMUX), on-line integration (INT), 
correlation computation (CORR) and on-line analysis (AN), as well as data 
telemetry (TEL) and remote control (CANA). The latter additions of remote 
control and data transfer are usually needed when such a radar system would 
be used in operational applications as a wind profiler. 


ANTENNA 



Figure 25. Transmit- receive duplexer (also called T/R switch) 
using 3-dB hybrids and pin-diode switching. 
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receiver hybrid and exits through the resistor port. In the reception case the 
diodes are open and the split power from the antenna is combined In-phase In the 
receiver hybrid to pass through the receiver port to the receiver Input. Any 
noise, which is still generated in the transmitter during the reception phase is 
combined Into the resistor port of the receiver hybrid. The response and re- 
covery time of the duplexer should be In the order of the range sampling time, 
l.e. typically less than 10 ps. To insulate a peak transmitter pulse amplitude 
of 10 5 W to a fraction of Watt, which will not affect the receiver, the decoup- 
ling attenuation has to be better than 60 d&. These specified values can be 
obtained with the described pin-diode hybrid switch. Other kinds of transmit- 
recelve duplexers, using quarter-wave coax lines or parallel and serial reso- 
nance circuits are also used. These are not so suitable because of the resonance 
circuits, which cause ringing effects and affect the data in the first sample 
range gates. If circular polarization is applied (in IS radar operations) the 
duplexer can be replaced by a hybrid, which provides some 30 db decoupling 
between the transmitter and receiver port. The remaining RF signal leaking 
through the receiver port is usually attenuated by a so-called receiver protec- 
tor, which is a fast pin-diode switch disconnecting the receiver input from the 
hybrid receiver port during the transmission phase. In order to further protect 
the receiver during eventual malfunctions of the duplexer or the receiver pro- 
tector the status of these devices is monitored by proper hand-shake and inter- 
lock systems, which prevent the transmitter to be turned on in case of a fai- 
lure. 


ARRAY] 



Figure 26. Block diagram of the Chung-Li VHF radar, which bases 
on the principal lay-out of Fig. 24 (from Brosnahan et al . , 
1983). 
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In Fig. 26 a block diagram of the Chung-Li VHF radar Is shown, which is a 
dual-mode radar. Dual-mode means that the radar can be used In the Doppler as 
well as in the spaced antenna mode (for velocity and Interferometer measure- 
ments). For the spaced antenna mode It has three separate antenna arrays for 
reception which are combined in phase for transmission. In the Doppler mode the 
antenna arrays can be either combined to point into the same direction during 
transmission and reception, or they can point at the same time into three dif- 
ferent directions. The beam steering is done by phase shifters and the power 
distribution by power dividers, which are made from coaxial cables. The trans- 
mlt-receive duplexer is called T/R switch (transmit- receive switch) in this 
diagram. Three transmitters are feeding the three antenna arrays. For providing 
the essential phase equality at the three antenna arrays, the transmitter drive 
RF signals are phase-controlled. In order to provide the phase-coherent radar 
system, the same oscillator controls the transmitters as well as the three 
receiver channels. Following the sample- and- ho Id (S/H) circuit and analog-digi- 
tal converters (A/D) the data are preprocessed, i.e. coherently integrated, and 
transferred to the mini -computer. From there the data are dumped on tape. By the 
computer also the radar operation Is controlled, particularly by loading the 
system synchronizer (radar controller) and starting the operation. The system 
synchronizer also steers the antenna beam directions. We now will briefly de- 
scribe in more detail some typical sub-units of such a radar system. 

3.1. Transmitters, Receivers and Oscillators 

In Fig. 27 the principle block diagram of a transmitter is shown. The os- 
cillator generates a continuous wave signal, which is pulse modulated (in ampli- 
tude and/or phase) by a fast switch controlled by pulse trains from the radar 
controller. The switching needs to be done with transition times, which are 
reasonably faster than the duration of the shortest pulse in single pulse appli- 
cations or the pulse element in phase coding. Thus the switching of 1-us pulse 
elements needs to be done within 100 ns or so. The pulsed RF signal, usually at 
a level of several mW, is amplified in driver stages up to peak power levels of 
several kW. The final amplifier of NST radars usually yields several ten kW peak 
RF output power. Most frequently vacuum tube amplifiers are applied in the high 
power stages, which need power supplies with at least 5-10 kV high voltage. 
Semi-conductor devices needing much lower voltages are now getting in more use, 


T R A N S H I T T E R 



Figure 27. Principle block diagram of a typical MST 
or incoherent scatter radar transmitter. 
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particularly when distributed antenna modules are employed Instead of one large 
antenna array. Because of the high output power requirements of more than 1 MW, 
the incoherent scatter radars apply klystron tubes in the final transmitter 
stage, which need pulsed power supplies with 50-100 kV high voltage. The high 
power - high voltage power supplies usually are equipped with a very fast pro- 
tection switch (a so-called "crow bar"), which discharges the power supply 
within small fractions of a second in case of a transmitter malfunction to 
protect the final power amplifier tubes from damage. 

Following the final amplifier the transmit-receive duplexer switches the 
antenna either to the transmitter or to the receiver as described in the last 
chapter. To allow a measurement of the output power a directional coupler is 
used after the duplexer output. By measuring the forward as well as the re- 
flected power with two directional couplers the matching (i.e., the voltage- 
stand ing-wave- ratio: VSWRJ of the antenna is usually monitored. If a phase 
controller is needed to provide exact phasing of separate antenna arrays, it 
usually is inserted in the low power stages of the transmitter. Such a device is 
principally described by the block diagram of Fig. 28. The signal from the 
transmitter output, measured with the directional coupler, is compared with the 
local oscillator CLO) reference signal with specified phase (for multnple 
transmitter beam steering). The phase offset between both these signals is 
sampled and used to control the phase shifter, which changes the phase of the 
input signal (RF*r> ). The RFout signal, available with the corrected phase at 
the output of the phase controller, Is then fed to the next amplifier stage in 
the transmitter. 

Usually the local oscillators are more complex because of necessary fea- 
tures, which are outlined in Figs. 29 and 30. For a super-heterodyne receiver, 
i.e. a multiple mixing of the receiver input signal to one or two intermediate 
frequencies with subsequent mixing to baseband, a main oscillator in any in- 
stance allows for the phase synchronization between these different frequencies 
(Fig. 29). In case of multi-frequency channel operation, also a set of phase- 
locked frequencies has to be generated. In practice the transmitter as well as 
the receiver oscillator signals have to be further processed as shown in Fig. 
30. The main local oscillator signal is divided into two signals in the trans- 
mitter and the receiver path, respectively. In the transmitter path frequently a 
phase flip between 0° and 180° is implemented to allow phase coding (phase 
modulation), followed by the switch to turn on and off the transmitter signal 
(pulse amplitude modulation). In the receiver path a switch turns off the local 
oscillator signal during transmission to reduce any feed-through of leakage of 
the transmitter signal into the receiver. After power dividing of the local 
oscillator signal one path is phase shifted by 90° to allow baseband mixing to 
the quadrature components, i.e. real (Re) and imaginary (Im) part of the base- 
band signal. The control signals ( FLP = flip, TXP = transmitter pulse, and LOP = 
local oscillator protect) are generated by the radar controller (see Chapter 5). 


Figure 28. Block diagram of 
an automatic transmitter 
phase-controller, which is 
needed to phase-synchromze 
multiple transmitters feed- 
ing parts of an antenna. 
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Figure 29. Slock diagram of an oscillator circuit for super-heterodyne 
receivers with two local oscillators LOl and L02 and the phase- 
coherent transmitter oscillator control. 
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Figure 30. Block diagram of the master oscillator to control the 
transmitter and the receiver, showing the particular phase 
shifting and switching circuits for the transmitter and receiver 
osci 1 lator signals. 
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Figure 31. Block diagram of the principle 
receiver circuitry of MST and IS radars. 

The receiver Is principally described by the block diagram of Fig. 31. The 
low-level signal from the antenna (either directly in case of separated trans- 
mitter and receiver antennas or via the transmit- receive duplexer in case of a 
common receiving-transmitting antenna) is fed to the low-noise preamplifier 
through a band-pass. Frequently also a calibrated noise Injection, which is 
turned on and off by the radar controller, is directly fed into the receiver 
front end. This is necessary to allow a well calibrated estimate of the receiver 
amplification factor and thus a deduction of the system temperature and the 
absolute power of the received signals. In the first mixer the receiver signal 
Is mixed with the signal from the first local gscillator (L01) to the inter- 
mediate frequency (IF). Again a pass-band filter separates the wanted IF from 
the unwanted mixing products. Before further amplification usually a variable 
attenuator is implemented to adjust the total amplification of the receiver and 
thus the output levels for the final baseband mixer and the analog-digital con- 
verters. Following the mixing to baseband with the L02 signal, i.e. the qua- 
drature detection (Re and Im), a lowpass filter Is applied which usually should 
match the signal bandwidth (given by the bandwidth of the transmitted pulse for 
MST radars and by the bandwidth of the scattering process for IS radars). An 
example of a low-noise receiver front-end of a UHF incoherent scatter radar, 
which allows the calibrated noise injection, the reception of circular and 
elliptical polarization as well as sideband conversion, is shown in Fig. 32. The 
polarization adaption is achieved by phase changes of the horizontal (H) and the 
vertical (V) signal components in the polarizer and the sideband conversion by 
mixing from the upper side (1053.5 MHz) or the lower side (813.5 MHz). Depending 
on the output from the last hybrid of the polarizer also the sense of rotation 
of the polarization (left- or right-hand rotation) can be selected. 
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120 1 15 MHz 


Figure 32. Block diagram of the 933-MHz receiver front-end of the 
EISCAT incoherent scatter radar receiving stations. The hori- 
zontal (H) and vertical (V) components of the signal from the 
feed horn of the dish antenna are combined by a polarizer (phase 
shift and hybrids) and mixed to the intermediate frequency 120 MHz. 
A noise generator, which is switched on and off by the radar con- 
troller, injects a calibration signal into the receiver front-ends. 


3.2. Examples of Radar Systems 

There are some technical system specifications, which principally apply for 
all atmospheric/ionospheric radar systems. A block diagram of a typical trans- 
mitter/receiver system of an MST radar is shown in Fig. 33. It depicts the 
design of the portable SOUSY Mini-VHF-Radar (46.8 MHz), which was operated at 
the Arecibo Observatory in 1980. According to individual requirements, many 
variations and modifications of such a system can be done, but Fig. 33 describes 
a generally basic design. The 120-MHz signal of the main or paster oscillator 
(MO) is divided by four to obtain a 30-MHz Intermediate frequency (IF) which can 
be modulated (MD) in amplitude (RFC = R£ control) and phase (FLP.= flip). This 
is to provide phase coding and DC-elimination (see Chapter 4 for more details). 
A similar divider generates the 0° and 90° signals, which are necessary for qua- 
drature detection. The operational frequency 46.8 MHz is generated by mining (X) 
with the local oscillator (LO) signal at 76.8 MHz. It is amplified in the trans- 
mitter (TX) and fed through the transmit-recelve duplexer (TRX) to the antenna. 
The signal received at 46.8 MHz is gated (RGT = receiver gate), amplified In the 
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SOUSY MINI VHF-RADAR MPAe -LINDAU 



TRX RFC FLP Re Im RGT 


Figure 33. Block diagram of the transmitter-receiver system of 
the portable Mini SOUSY-VHF-Radar , which was used with the 
Arecibo Observatory dish antenna for MST studies. 

receiver (RX), mixed with the LO-signal to an intermediate frequency (IF) of 30 
MHz, and mixed down (X) to the baseband in the quadrature detector. The two 
quadrature outputs, the .real part (Re) and the Imaginary part (lm), are low pass 
filtered CLP) to match the receiver bandwidth to the bandwidth of the trans- 
mitted pulse. 

The transmitter peak power typically is between 1 kw and 1 MW for MST 
radars with duty cycles up to several percent, while incoherent scatter radars 
use duty cycles of some 10% and peak powers above 1 MW. The transmitter band- 
width must cover the shortest pulse length of 1 us (*2 MHz RF bandwidth). The 
transmitter is normally operated in class-C mode to achieve an optimum effi- 
ciency. Incoherent scatter radar transmitters use mostly high power klystrons. 

The receiver linearity range usually exceeds 60 - 80 dB in order to avoid 
saturation with strong clutter signals or in the case of large dynamic range 
variations of the atmospheric signal. The phase flip and the quadrature detec- 
tion must be within less than a few degrees accuracy, and the amplitude ratio of 
the quadrature components must not deviate from unity by a few percent in order 
not to distort the Doppler spectrum analysis. The stability and the phase noise 
of the oscillators is usually better than 10" 4 Hz/ms to allow a good accuracy of 
coherent detection. The receiver base-bandwidth (postdetection) has to be about 
1 MHz to provide the resolution of the shortest pulses of 1 us. The intermediate 
frequency bandwidth needs to be at least twice as large as the base-bandwidth. 
Bessel filters are usually applied after the quadrature detector as low-pass or 
post-detection filter in MST radars, but Butterworth filters in IS radars. 

The receiver noise figure of an MST VHF radar need not be better than a few 
dB (some 100 K), since the sky noise level (> 1000 K) determines the sensiti- 

vity. A noise calibration signal of 1000 K, say, should be injected into the re- 
ceiver frontend in order to provide an absolute power calibration (which also 
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Figure 34. Sky noise temperature as function of frequency 
and antenna elevation angle (from KRAUS, 1966). 



Figure 35. Sky noise map measured as function of declination and 
right ascension on 64 MHz (after KO, 1958). 
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needs a continuous monitoring of the transmitter power). However at higher 
frequencies in the UHF band the sky noise is no more the determining factor as 
can be seen from Fig. 34. Therefore above several 100 MHz (e.g. in case of the 
Incoherent scatter radars or for the UHF ST radars where the signal-to-noise 
ratio is usually very low) the receiver frontend has to be carefully designed 
for lowest equivalent noise temperatures. Usually system temperatures of several 
ten Kelvin can be achieved at some incoherent scatter UHF radars. The sky noise 
(from galactic radio sources) changes as function of time of the day and antenna 
pointing angle as depicted by Figs. 35 and 36 for frequencies of 64 MHz and 
53.5 MHz, respectively. This changes the signal-to-noise ratio even if the 
reflectivity of the radar volumes would be constant. The known variation of sky 
noise and in particular the well known ephemeris and flux density of point-like 
radio sources (e.g. Cassiopeia or Cygnus) yield on the other hand a well suited 
means to calibrate antenna diagrams and antenna gains. 



Figure 36. The diurnal variation of sky noise level on 
30. /3i. August 1978 measured at the antenna of the 
53.5-MHz SOUSY-VHF-Radar (52°N) in West Germany. 

The equivalent noise figure of 0 dB corresponds to 
a system temperature of 290 Kelvin. 


In Figs. 37, 38 and 39 we show more block diagrams of three other radars, 
the Poker Flat MST Radar operating on 50 MHz, the MU radar operating on 46.5 
MHz, which is used for MST as well as incoherent scatter observations, and the 
EISCAT radar systems, which operate as incoherent scatter radars in the 933 MHz 
and the 224 MHz bands. In Fig. 40 a block diagram of a 915 MHz wind profiler 
radar and in Fig. 41 the block diagram of the Arecibo bistatic S-band (2380 MHz) 
CW-radar Is shown. The latter radar operates on 2380 MHz in phase-coded continu- 
ous wave (CW) mode. The experienced reader will notice several of the basic 
features in these schematic diagrams which we have mostly outlined in the fore- 
going chapter. More technical details of these radars should be found in the 
relevant literature. 




COAXIAL ANTENNA ARRAY 


Figure 37. Block diagram of the Poker Flat 
MST Raaar (from BALSLEY et al., 1960). 
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Figure 38. block diagram of the MU radar 
in japan itrom FUKAO et al., 1985a, b). 
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Figure 40. Block diagram of a 915-MHz Doppler radar 
wind profiler (from FRISCH et al., 1986). 
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Figure 41. Block diagram of the bi-static S-band (2380 MHz), 
phase-coded CW (continuous wave) radar system used at the 
Arecibo Observatory in Puerto Rico for high resolution 
stratospheric turbulence studies (from WOODMAN, 1980). 


4. RADAR SIGNAL ACQUISITION AND PREPROCESSING 

The simplified schematics of radar systems, shown in Figs. 6a and 6b, are 
used for a basic explanation of the radar operation and data-acquisition proce- 
dure. After we have briefly explained the lay-out of antennas, transmitters and 
receivers, in the last chapters we will now explain the formalism of a basic 
coherent system, the quadrature detection, the digital sampling as well as 
subsequent data acquisition and phase decoding procedures. For more details see 
WOODMAN and GUILLEN (1974), ROTTGER and SCHMIDT (1979), SCHMIDT et al. f (1979), 
CARTER et al . , (1980), CLARK and CARTER (1980), SATO and WOODMAN (1982), RASTOGI 
(1983), for instance. Since these signal processing steps are mostly done in 
specially designed hardware processors, we now briefly summarize the underlying 
concepts. 

4.1. Coherent Detection 

The oscillator generates a signal s° at the angular frequency Qo = 2nfo 
where fo = c/0o is the center radar operation frequency. A pulse tram, gener- 
ated by the radar controller, imposes a modulation to this signal. After am- 
plification in the transmitter (TX in Figs. 6a and 6b) the radar signal 

s(t) = a(t) exp(i(Qot + Pit))) (5) 

Is transmitted, where a(t) determines an amplitude modulation (by the pulse 
train) and 0(t) corresponds to a phase modulation (for coding), and i = v-1. 
a(t) and p(t) are slowly varying as compared to Qot. The radar signal is scat- 
tered/reflected from the radar volume and reaches the receiver (RX) via the same 
or a separate antenna. Additionally, noise (sky noise and interference) Is 
received and adds to the radar echo. 


Q'ZL 
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For simplifying the explanation we will assume for a while that the trans- 
mitted signal s(t) Is not modulated, l.e. a(t)=const and 0(t)=const, and the 
echo results only from a small radar volume at a given range. The echo signal 
s’(t) plus noise r(t), which are both band-limited because of the scattering 
process and/or the receiver bandwidth, can be represented by 

c’(t) = s’(t) + r(t) = ai(t)cos Qot + ia 2 (t)sin Qot, 

where ai(t) and a 2 (t) are independent Gaussian variables in a pure scattering 
process. The uncorrelated noise r(t) contributes only uncertainties to these 
estimates. After linear amplification In the receiver, c'(t) is coherently 
detected by multiplicative mixing with s°. After low-pass or post-detection 
filtering (to eliminate high frequency components 2Qo, which are generated 
during the mixing procedure), this yields 

c ( t ) = a + (t)cos 0 + (t) + iaMtlsin 9Mt), 

where a + (t) = (ai 2 (t) + a 2 z(t))^V2, (6a) 
e*(t) = arctan (a 2 (t)/ai (t ) ) . 

The phase 3+(t) = Qot - 9’(t) is given by the Doppler frequency Qd = -4nWfio, 
which is due to the bulk motion V’ of the scatterers. The time variable phase 
e’(t) = e ( t ) + d(t), where 8(t) is caused by the fluctuations of the scatterers- 
/ref lectors In the radar volume. The amplitude a + (t) is a measure of the reflec- 
tivity of the scattering/reflection process. These latter statements are only 
valid if the noise contributions are separated from the signal. The coherently 
detected complex signal (+ noise) can be expressed In the form 

c(t) = x ( t ) + iy(t) T (6b) 

where the real part x(t) = a + (t) cos 8 + (t) is called the m-phase component, and 
the imaginary part y(t) = a + (t) sin 9 + (t) is called the quadrature component. 
Both components, x and y, are called the quadrature components. The Fourier 
transform of c(t) is 

X(Q) = |c(t)exp(-iQt)dt = 5KQ) + i?lQ), 

which yields the periodogram P(Q) = X 2 + 9 2 i it is often also called power or 
Doppler spectrum. The measured P(Q) is the convolution of the spectrum of the 
refractive index or the reflectivity fluctuations in the radar volume with the 
spectrum Pt(Q) of the transmitted wave form (for the common case of pulsed 
transmissions), multiplied by the bandpass characteristics Pr ( Q ) of the recei- 
ver. Since P(Q) is much narrower than the envelope of Pt(Q) and Pr(Q) in MST 
radar investigations, these instrumental effects can mostly be disregarded. 
However, in incoherent scatter applications the spectrum of the scatter process 
P R (Q) is mostly wider than the envelope of Pt(Q), which leads to the fact that 
the instrumental effects of receiver and transmitter pulse bandwidth have to be 
accounted for during the signal analysis. 

4.2. Digital Sampling 

We now take into account that the signal c(t) results from scattering and 
reflection in certain altitude ranges. Frequently these signals which originate 
from scatter or reflection of a transmitted radar signal are called "radar 
echoes". We now have to note that the transmitted signal is modulated in order 
to resolve the range from where the radar echoes arrive from. As sketched in 
Figs. 1 and 2, the echoes from different ranges then occur at different times at 
the analog receiver output. In order to allow adapted signal processing of 
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Figure 42. Range-time-amp 1 i tude diagram of single-pulse 
transmitter modulation and coherent MST radar signals. 

digitized data, the signal c(t) is sampled at discrete time intervals tk = k6t« 
In the §nalog-£1gital converter (ADC). This is illustrated In Fig. 42, which 
shows the amplitude variations of one of the quadrature components of a coherent 
MST radar signal as they could be monitored with an oscilloscope connected to 
the output of the postdetection filters, which follow the quadrature detectors. 
The sequences n = 1,2,-.- can be assumed to represent successive oscilloscope 
beam deflections (along the vertical axis), triggered by the leading edge of the 
transmitter pulse. This pulse Is strongly attenuated by receiver gating, as was 
explained in the last chapter. It is flipped in phase by 180° from one pulse to 
the next (change from positive to negative amplitude), for the reason of DC- 
eliminatlon (see details in the following Chapter 4.3). Corresponding to the 
transmitter-pulse phase 0, the sign of the radar echo also changes from one to 
the next interpulse period. For convenience the phase of the radar echo is shown 
here to be similar to the transmitter-pulse phase. The echo phase can take any 
values, however, depending on the length of the phase path from the transmitter 
to the radar volume and depending on the relative changes of the scatterers 
within the radar volume. 

On the vertical axis of Fig. 42 the time tk is given, which directly can be 
converted to range by means of the definitions of Chapter 1.1. The sampling time 
Interval 5t» should be equal to the radar transmitter pulse length 5tt , since 
this usually (in case of volume filling scatter) yields an optimum matching to 
the range gate width or resolution 6r. However, in certain cases when narrow 
reflecting layers shall be resolved, shorter sampling intervals 6t* < 5t r can Pe 
useful; this Is called "range oversampling". 

The signal and the noise had passed the receiver and postdetection filters 
and are therefore bandlimited. The response time of the receiver, which is 
approximately inversely proportional to the filter bandwidth, should also be 
equal to 6t», respectively 6tt . The subscript k of tk is the serial number of 
the range gate with K = 0,1,..., K-1, where k = 0 corresponds to the beginning 
or a certain level of the transmitter pulse. Since the quadrature components are 
digitally sampled, we can write equation (6b) in the form Ck = Xk + lyk. We also 
call Ck the complex raw data samples. 
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The number of sampling time steps between successive radar transmitter 
pulses Is K, which is also the number of sampled range gates. The interpulse 
period Is Tipp = K 6U. The cycle k = 0, K-1 is repeated once with 

k’ = 0, K’-1 with a reversed phase of the transmitter pulse. Both cycles, 

denoted by k and k’ (with k’ 2 k and K’ 2 K), determine one radar cycle 2*Ti PP 
(note that here the radar cycle is two times the interpulse period, because of 
the introduced phase flip). The serial number of radar cycles is given by 
n = i , . . . , n. One radar burst is determined by N radar cycles, which last for ti 
= 2N Tipp . The generation of all pulse trains, needed to control these cycles of 
the transmitter and receiver-ADC-integrator system, is done in the radar con- 
troller, as described in Chapter 5. 

4.3. Coherent Integration and Preprocessing 

We will discuss here a standard preprocessing procedure of the pulse-to- 
pulse technique applied with MST radars. The spectrum which one would obtain for 
a coherent echo (HST radar) with the sampling rate Ti PP is very wide and mostly 
contains high frequency noise power. The signal power Is confined to relatively 
low frequencies only (f = 1/t€ << 1/2T< PP ), where xe is the signal coherence 

time. It is evident therefore that low pass filtering of a coherent signal, done 
before the spectrum analysis, will not change the signal characteristics but 
eliminates high frequency noise contributions. The simplest form of low pass 
filtering is just the complex addition of the signal + noise samples over an 
interval u << xc . A readily noticed effect of such an averaging over N inter- 
pulse periods is the reduction of the number of total samples by a factor 1/N. 
We have to note that this processing step of coherent integration cannot be 
applied for Incoherent signal processing, because the signal coherence time is 
shorter than the interpulse period. This also means that the described DC- and 
clutter-elimination is not applicable for incoherent signals. The standard 
procedure of autocorrelation function and spectrum computation, however, is 
compatible in both applications. 

Since the noise (r) and the signa! (s) are independent of each other, their 
quadrature components add to Ckn = Ckn r + Ckn*, where Ckn r ( = Xkn r + iykn r ) 
are the quadrature components of the noise and Ckn* those of the signal. A DC- 
bias can result from a constant, instrumental ly introduced voltage at the recei- 
ver output Ckn 1 , or due to radar clutter Ckn c (clutter = echo from a fixed 
target). All these contributions are additive: 

Ckn = Ckn r + Ckn 8 + Ckn 1 + Ckn c . 

We have to take Into account that each radar cycle yields two samples per range 
gate, namely k and k\ The samples of signal and clutter are shifted by 180° 
from k to k 1 , because the phase 0(t) of the transmitter was flipped by 180° 
(change 0(t) by 180° in equations (5) and (8)). This can be accounted for by 
changing the sign of Ck-n when averaging 


Ck , k ’ = - I (Ckn “ Ck ' n ) 

N , 
n= 1 

N 

= - 2 (Ckn r _ Ck’n r + Ckn® “ Ck’n® + Ckn 1 - Ck ’ n i + Cknc “ Ck'ncJ. 
n- 1 

Since for the Instrumental DC-bias Ckn 1 = Ck’n 1 , it is eliminated by averaging. 
This Is called Instrumental-DC elimination (DC = direct current, better to say: 
constant voltage contribution). 
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Because of the transmitter phase flip, we find that: Ckn* = - Ck’n* and 
Ckn c = - Ck»n c . Since the noise Is independent from one to the next interpulse 
period, a change in sign of Ckn r does not change its statistical properties. We 
thus obtain, since k* = k: 


2 N 

Ck = - I (Ckn r + Ckn 1 

N , 
n- 1 


+ Ckn c ) 


(7} 


This averaging, commonly called "preintegration" or better "coherent integra- 
tion" , has become a standard process in MST radar operations. It yields cohe- 
rently integrated data samples Cki (1 = 1 L) at the time 

ti = 2 - 1 * N - K - 6t* 

for the coherent integration period ti = 2-NK-6U = 2-NTi pp . 


If the summation in equation (7) extends over a time period ti = 2N-Ti PP , 
which is much longer than 6t« , the high-frequency noise contribution vanishes. 
For a coherent integration period tL = 2-L*NTi PPt with (L >>N), which is much 
longer than xe , the signal contribution approaches zero since it is fading in 
amplitude and phase. Only the clutter contribution 

Ck c = j- I Cki 
L 1 = 1 


remains, since it is constant in amplitude and phase. This can be used to elimi- 
nate the clutter component by means of a digital high pass filter operation: cki 
= Ck i - Ck c , This operation Is called clutter-OC elimination, and is done after 
the coherent integration. 

For the coherent integration given by equation (7), the number N of added 
samples has to be selected carefully. It is evident that the integration period 
has to be much shorter than the typical time scale of signal variations due to 
the fluctuating scatter process as well as due to the frequency changes resul- 
ting from the bulk motions of the scattering or reflecting medium. The advantage 
Is, however, that the processing can be made very efficient since only additions 
(subtractions) but no multiplication of the raw data series with a weighting 
function are necessary. It is also required that the real part x and the imagi- 
nary part y of the detected signal are correctly in quadrature (orthogonal). 
Ideally they must be phase-shifted by exactly 90° and must have equal amplitudes 
on the average, otherwise a distortion of the Doppler spectrum results. Accu- 
racies of less than a few degrees phase difference and less than a few percent 
amplitude difference are tolerable, however, and can be obtained by proper hard- 
ware adjustment. 

The coherent integration of the quadrature components, formulated by equa- 
tion (7), Is normally done in a digital preprocessor, called adder or Integra- 
tor. Since this coherent integration is a low pass filter process it can be done 
also in an analogue filter, such as the clutter elimination which is a notch or 
high pass filter operation. Obviously the digital processing Is much more versa- 
tile and flexible. Both analogue and digital preintegration reduce the number of 
data samples by some orders of magnitude, and compress the huge data flow from 
the ADCs to make it manageable for the host computer. This Is the evident advan- 
tage of this preintegration process. We have to note, however, that this cohe- 
rent integration process looses its advantages for higher frequency 0500-1000 
MHz) radar applications. The reason Is that the signal coherence time Is In- 
versely proportional to the radar operating frequency and consequently the 
coherent integration time has to be reduced. It is often thought that this 
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process of coherent integration also Increases the signal-to-noise ratio since 
the voltages of the coherent signal but the power of the incoherent noise add. 
This leads to an Improvement of the signal-to-noise (power) ratio by the factor 
N, since the noise bandwidth is reduced by the factor 1/N, whereas the signal 
bandwidth remains unchanged. If one reasonably defines the noise in such a way 
that its bandwidth Is equal to the signal bandwidth (e.g., ROTTGER and LARSEN, 
1989), the coherent integration process only reduces the wide noise bandwidth to 
the acceptable limit close to the signal bandwidth, and the signal-to-noise 
ratio does not change by this coherent integration. The advantage of the coher- 
ent Integration process Is still the essential reduction of the number of raw 
data samples, without giving away information on the signal. 



Figure 43. Range-time-amplitude diagram of a long-pulse 
IS radar modulation and Incoherent signals. 

This kind of coherent integration used in the pulse- to- pulse technique 
cannot be applied for incoherent scatter signals, because the signal coherence 
time Is shorter than the interpulse period, and thus the preprocessing needs to 
be done differently as was principally explained already in Fig. 3. For clear- 
ness the same example of an Incoherent signal (Fig. 43) is drawn In the same 
manner as for the coherent MST radar signal (Fig. 42). We notice in Fig. 43 the 
"Incoherence" from one to the next interpulse period, which leads to the fact 
that the correlation-function and spectrum analysis has to be done within one 
Interpulse period and the ACFs and spectra then are Integrated subsequently. The 
Fig. 43 shows that a fairly long pulse is still needed to obtain samples for the 
longest lag of the autocorrelation function (ACF) In order to achieve the neces- 
sary frequency resolution. This obviously deteriorates the range/altitude re- 
solution. To overcome this problem, so-called multi-pulse modulation Is applied, 
and is explained In Fig. 44. This modulation is most efficiently applied by 
modulating different frequency channels F1-F6 (within the passband of the trans- 
mitter and receiver) in order to fill the gaps between the sub-pulses of this 
pulse pattern. This Increases the average power to the maximum allowable level, 
but needs a particular multi-frequency channel design of transmitter and recei- 
ver. To Improve the range resolution further the subpulses can be Barker-coded. 
Such a multi-pulse multi-frequency Barker-code scheme is successfully used in 
the El SCAT Incoherent scatter radar systems (e.g.. TURUNEN, 1986). 
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Figure 44. Multi-pulse multi-frequency-channel pulse pattern used 
In incoherent scatter applications to improve the range reso- 
lution and allow measurements of the lag products of the signal 
autocorrelation function for the temporal lag Indices j = 1 to 6 
(after TURUNEN, 1986). 

The application of this modulation scheme as well as the successlng "alter- 
nating code scheme" has become standard praxis in incoherent scatter work as has 
the "complementary code scheme" in the MST radar work. We will explain the 
latter coding scheme after we have briefly outlined some basic formulations of 
the correlation analysis, which Is applied In principle also in the decoding 
procedures. Again we have to consider that most of the described preprocessing 
can be done by specially designed digital hardware. 

4.4. Correlation and Spectrum Analysis 

After the preprocessing by appropriate coherent integration (for coherent 
signals) or immediate correlation function processing (for Incoherent signals) 
one usually deducts specific signal parameters, such as signal power, Doppler 
shift and spectrum width. Usually these are the standard parameters in MST 
applications, whereas in incoherent scatter applications more information Is 
extracted from the spectral shape (see BARON, 1977) by special fitting routines. 
Before doing this, one has to eliminate the noise which is remaining within the 
bandwidth given by the signal. The elimination of the noise and further signal 
processing can be done in the time domain by correlation (covariance) function 
analysis or in the frequency domain by spectrum analysis. These data processing 
steps, which follow the hardware preprocessing, can be suitably done by speci- 
ally designed software in real time In. the host computer or off line In separate 
computers by reading the data from tape. In incoherent scatter applications 
usually the correlation function computations are done in a special purpose 
hardware correlator, whereas In MST radar applications the essential preprocess- 
ing procedures, which need special hardware, are the coherent integration and 
decoding. In the latter case further processing to obtain the ACFs or spectra 
need to be done after the preprocessing. 

Since the covariance function and the power spectrum are Fourier transforms 
of each other, both contain the same relevant information. Depending on the 
purpose and the feasibility of the analysis, either of both approaches, Is used 
In practice (RASTOGI and WOODMAN, 1974; WOODMAN and GUILLEN, 1974; Hagfors, 
1977; ROTTGER and SCHMIDT, 1979; SCHMIDT et al., 1979; CARTER et al., 1980; 
CLARK and CARTER, 1980; SATO and WOODMAN, 1980; RASTOGI, 1983). 
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The complex autovariance function of the quadrature components c(t) Is 


?(t) 


c(t) -c*(t + x)dt. 


where x Is the temporal displacement, and the * denotes the complex conjugate. 
In digital form 


1 

Rk(ij) = “ I Cki'C*ki+j = CkTC*kHj, with 3=0 J-1 and J < L, 

1 j 1 = 1 

where j Is the lag Index defining the lag Xj . For a fixed range k, we find ci = 
xi + lyi, and the covariance function becomes 

R(xj) = (xixi+j+yiyi4j) + 1 (xi + jyi-xi yi+j ) 

= Rr (Xj ) + i Ri ( Xj ) , 
or R(xj ) = |R(Xj ) | exp(1(Xj)), 

with | R | = (Rr 2 +Ri 2 )i and 0 = arctan (Ri/Rr). 

In radar applications the term correlation function Is often used for R(Xj) 
instead of covariance function. The correct definition of the correlation func- 
tion g(xj) is given by the normalized covariance function: 

g(Xj ) = R(xj )/ R(0) . 

As well as the autocorrelation function of the series c(t) we also compute 
cross-correlation functions for two different series, ci(t) and C 2 (t), in the 
spaced antenna applications and use their modulus in the drift analysis and 
their phase and amplitude in the Interferometer analysis. In a principally 
similar way also cross spectra are computed from two time series. 

The power spectrum is the Fourier transform of c(T). weighted by W(T): 

P (Q ) = |* W(x ) g(x) exp(-iQx)dx 

or in digital form: 

1 J_1 

Pk(Qm) = - I Wj-Rkj exp(-iQ«jti ) , 

J 3=0 

where Wj is an arbitrary weighting function (e.g., Wj = 1 for all j, etc.) and 

Qm = m n/J • ti ; m = 0 J. The length of the Integration time interval is 

tn = 2NT*pp In the case of a coherent signal. For an incoherent signal, ti 
corresponds to the time interval over which lag samples of the ACF are collec- 
ted. 


4.5. Phase Coding 

Another preprocessing step which is generally similar to preintegration is 
the decoding procedure. The principal reason for applying coding/decoding (pulse 
compression) is to achieve a maximum average power at optimum resolution and 
maximum unambiguous range. A fairly low duty cycle of 0.1% and hence a low 
sensitivity would for instance result from using a single pulse of 1 ps duration 
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(6r = 150 m) and an interpulse period of 1000 us (r»*x = 150 km). Increasing the 
pulse length to 32 us would Increase the duty cycle to 3.2%, but deteriorate the 
range resolution to 4.8 km. The range resolution of 150 m can still be achieved 
by phase-coding the transmitter pulse In time (lag) Increments of, say 1 us. The 
decoding has to be done by cross correlating the received complex time series ck 
( separately for real and imaginary part) with the transmitted code s»: 

M-1 

Ck*' = I Ck*+« ■ Si (8) 

m=0 

where M Is the length of the code and k* corresponds to the range gate k at lag 
zero where the correlation function Ck*’ has a maximum. Thus, the decoding is 
nothing else but an integration of the signal amplitudes over several range gate 
samples, which are multiplied by the weighting factor s>. 

The simplest and most versatile phase coding scheme is the binary code, 
where the phase Is flipped between the two states -1(=0°) and +1(=180°). Since 
here the weights s* are +1 or -1, the multiplication needed in the decoding 
process reduces to an addition/subtraction operation. For this reason the decod- 
ing can be done in a preprocessor similar to a hardware preintegrator, or the 
integrator/decoder can even be one unit, as we will outline in Chapter 5. Note 
that the decoding has to be done for the amplitudes, since these contain the 
phase information An advantage In MST radar applications is also that the deco- 
ding/integration processes are Interchangeable, which reduces the number of 
operations by about two orders of magnitude (e.g., WOODMAN et al., 1980). 

The best codes for radar applications obviously are those where the range 
sidelobes of the correlation function (8) at k* if k are a minimum. Reasonable 
sidelobe suppression is achieved with Barker codes, where the phases are flipped 
in a sequence (e.g., +++++--++-+-+ for the 13-baud Barker code). The correlation 
function is 13,0,1,0,1,0,1,0.1,0,1,0,1. The best range sidelobe suppression, 
achievable with a Barker code, is M - 1 . The reason is that these sidelobes con- 
tain power from other range gates k* i k, which cause ambiguities. The sidelobes 
of the code extend out to ranges Mot* on both sides of the center peak. Since 
the codes cannot be infinitely long, because the minimum range is extended with 
the length of the code (which obviously Increases the minimum range from where 
echoes can be received), the sidelobe suppression of Barker codes is limited. A 
good sidelobe suppression is for instance needed if strong echoes occur only in 
a limited number of range gates. 

A better sidelobe suppression can be gained by application of quasi-random 
codes, which however need decoding before the pre integration as well as to 
transmit a long series of different code sequences. Also alternating codes 
appear useful in Incoherent scatter applications (e.g., LEHTINEN and HAGGSTROM, 
1987). The most suitable codes for MST radar work are the complementary codes 
(e.g., SCHMIDT et al., 1979), which theoretically have no sidelobes. A comple- 
mentary code consists of a pair of two code sequences si and S 2 . These have the 
property that their correlation function sidelobes are exactly equal, but oppo- 
site in sign. Normally both code sequences are transmitted at one and the next 
interpulse period, and the range samples of these two periods are preintegrated 
and decoded separately. The coherent addition of the decoded sequences then 
yields the total elimination of the sidelobes (k* =f k) if the signal is coherent 
from one to the next interpulse period. The zero lag value k* = k contains the 
total signal amplitude. As an example the 8-baud complementary code pair A, B: 

si - A = +++-++-+ and S 2 = B = +++ +- 



99 


yield the correlation functions: 

ci = 8,-1, 0, 3,0,1, 0,1 and C 2 = 8, 1 ,0,-3, 0,-1 ,0,-1 
and the addition yields: 

ci + C2 = 16,0,0,0,0,0,0,0. 

The peak value of this decoded function at k* = k is 2M, There is still the 
shortcoming of a long code, that it extends the shortest observable range. This 
can be overcome by transmitting a sequence of a complementary code and a short 
single pulse In one radar cycle. In Fig. 45 the range-time-ampl itude diagram for 
the simplest 2-baud complementary code with additional phase-flip for instrumen- 
tal DC-el imlnation is shown. All these decoding preprocessing steps are most 
suitably done In a special purpose hardware integrator-decoder. For coding 
purposes, the hardware radar controller needs to generate the necessary ampli- 
tude and phase control pulse trains for the transmitter. 



Figure 45. Range-time-amplitude diagram of a two-baud 
complementary code modulation with DC-phase flip 
used in HST radar applications. 


5. RADAR CONTROL AND DATA ACQUISITION 

We will here briefly describe the basic system units which are needed to 
process the digitized data in terms of coherent integration, correlation func- 
tion computation and decoding as well as the principle lay-out of a radar 
controller. 

We first describe in Fig. 46 the basic lay-out of an integrator or adder. 
The data from the ADC are range-gate by range-gate stored into memory. According 
to equation (7), the data from one interpulse periods are added onto the corre- 
sponding data from the successive interpulse period in the arithmetic-logic-unit 
(ALU). This is recurrently done as long as the integration period is chosen, 
i.e. the number of interpulse periods preselected for the coherent integration 
(controlled through the host computer). After the completion of the integration 
cycle, the coherently integrated data are dumped via the di rect-memory-access 
(DMA) to the host computer. 
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Figure 46, Block diagram of an 
integrator/adder preprocessor 
for coherent signal acquisition. 


The coherent integration can also be done in the same unit as is used for 
decoding and the correlation function computation is done. Both the correlation 
and the decoding operation consist of a multiplication and an addition cycle, 
where the former needs complex processing and the latter needs to be done for 
the quadrature components separately. In Table 3 we have outlined the principles 
of the computation of correlation function, decoding and coherent integration in 
a graphic form, which should allow better understanding of the hardware configu- 
ration sketched in the diagram of Fig. 47. 

TABLE 3 


COMPUTATION OF CORRELATION FUNCTION, DECODING, COHERENT INTEGRATION 


1=0 { 23 . . l-l L 
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The square blocks In Table 3 Indicate data samples at range or time gates. 
The samples In the first row are supposed to be the original data. The hatched 
part of the diagram indicates data samples which are accessible at a register 
through which the original data are shifted. The correlation function is com- 
puted according to the formula A in Table 3. To compute the zero lag, the range 
(time) samples of the first (unhatched) row with indices j=J and its replica In 
the last row j=0 are multiplied for each range or time label 1, and then all 
products are summed up over 1. To obtain the first lag, the last but one series, 
which is delayed by one sample with respect to the last row. Is multiplied and 
added In the same way. According to this algorithm all the following lags are 
computed. In praxis, the last lag is computed first, since the data sample at 
j=j-1 is available first after the original data series had been shifted by one 
in a shift register (see explanation of Fig. 47 on the next page). 

In a very similar way the decoding is performed (formula B in Table 3), 
since it comprises a cross-correlation function computation which can be per- 
formed by the same hardware setup. The only difference to the correlation func- 
tion computation Is, that Instead of the data time series the code is used in 
the delayed data series. In the example of Table 3 a 5-baud Barker code is 
shown. It is to be noted, that the first completely decoded data point is avail- 
able only after 5 steps of j at J-5, when the 5 bauds of the code are completely 
existent in the shift register (for longer codes of course the corresponding 
longer delays have to be considered). The combination (1.e. # multiplication and 
addition according to the formula B) of series j=J-5 with j=J yields the decoded 
value at 1=0. As another example for instance j=0 with j=J yields the decoded 
sample at 1=3. Note that a further shift than given by 1=0 does not allow the 


CORRELATOR - DECODER - INTEGRATOR 



Figure 47. Block diagram of a hybrid version 
of a correlator-decoder-integrator. 
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complete decoding unless the data time series j=J Is further extended beyond 8 
samples in this example. In praxis this means that always a data set has to be 
available which is longer than a single data series by number Lc , where L c Is 
the number of bauds in a code minus one. Decoding of complementary codes is done 
In principally the same way, but as a final step the two cross-correlation 
functions resulting from the decoding of code si and Its complement S 2 have to 
be added. The coherent integration Is simply done, according to formula C in 
Table 3, In the same hardware processor by just adding the samples for each 1 
(or k) separately over a preselected number of interpulse periods. 

We have to note that all these procedures have to be done separately for 
the real part as well as for the imaginary part of the quadrature components. 
For convenience we explain only the integration, decoding and correlation func- 
tion computation for a single component. 

In Fig. 47 we show a block diagram of such a described hybrid "correlator- 
decoder-integrator". We admit that other realizations may be more practicable, 
but this layout should show that such a hybrid processor Is feasible In prin- 
ciple. The different modes can be understood by following the paths which are 
selected by different multiplexer’ positions. In the Integrator mode the multi- 
plexer A (MUX A) is connecting ports 1 and 3, multiplexer B ports 2 and 3, 
multiplexer 0 ports 1 and 3 and multiplexers D1 and D2 provide that each multi- 
plier Is connected to the correct range gate. The data are transferred from the 
buffer memory to the multipliers, where the sign is changed according to the 
corresponding phase flip of the transmitter pulse. The multiplier of course 
could be eliminated if provision would be made to allow also subtractions in the 
adder. For each range gate the data are added in the accumulator, from where 
they are multiplexed into the result memory after the specified number of cohe- 
rent integrations. Depending on the size of the result memory several sets of 
coherently integrated data are stored and afterwards dumped to the main compu- 
ter. 


For correlation function computations the multiplexers have to be in the 
positions A: 1-3, B: 1-3 and C: 1-3. The main part of the correlator is the 
shift register by which the data series is shifted according to the lag indices 
of the decoding or correlation function. Each data sample of the shifted series 
Is then multiplied with the data sample of the non-shifted data series. For 
decoding the shifted data series is the code, and for auto-correlation computa- 
tion the shifted series is the same unshifted data series. After this multipli- 
cation the data are accumulated, i.e, summed up, and after completion of the 
correlation analysis dumped into the result memory from where they are transfer- 
red into the host computer. The correct selection of .the lag products and their 
redistribution to the result memory is again provided by multiplexers D1 and D2. 

For decoding of the original data series, the multiplexer A is In position 
2-3, multiplexer B in position 1-3, and multiplexer C in position 1-3. The 
original data are existent at the multipliers and according to the position of 
the code In the shift register, the decoding is performed by multiplication. If 
the multiplexer A is in position 2-3, B in 1-3 and C in 2-3 a coherently inte- 
grated sample set could be also decoded in the same processor and afterwards 
dumped to the computer. It is noted that also a correlation function computation 
of decoded as well as coherently Integrated data could be possible if the pre- 
processed data are fed back by multiplexer C in position 2-3 to the result 
memory. Although the described layout of a multi-purpose integrator-decoder 
correlator explains the three main processing procedures of MST and IS radars 
and the solution looks feasible, It may in praxis be more realistic to use three 
separate units for Integration, decoding and correlation function / spectrum 
analysis. 
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TABLE 4 


RADAR CONTROL 


PRP 

TRX 

TXP 

FLP 

FR1 

FR2 

RGT 


_TL 


I LTUl 


— I L 


— f 


cycle 1 


LOP -f 
RIA — 
ADC — 
SGN -J 
CAL — 
ANT J~ 
AU X — 


TL 


Sl 


cycle 2 

fl 

J 1 

_J 1 

TUT 



n 


PRP = prepulse (trigger) 

TRX = transmitter-receiver duplexer 
TXP = transmitter on 
FLP = phase flip (coding) 

FR1 = frequency I 
FR2 = frequency 2 
ANT = antenna control 


RGT r receiver gating on 

LOP = local oscillator protect on 

RIA = reset integrator address 

ADC = analog-digital conversion 

SGN = sign for integration 

CAL = calibration (noise injection) 

AUX = auxiliary 


All the timing control of the outlined procedure results from the radar 
controller or system synchronizer, which of course also generates the control 
pulses for the transmitter, the receiver and the ADC. A typical example of pulse 
trains of control pulses, which are generated by the digital radar controller, 
is demonstrated in the example of Table 4. Following a pre-pulse (PRP), which Is 
used to to trigger external control or monitor devices, the transmit-receive 
duplexer (TRX) is switched on, followed by a receiver gate pulse (RGT) and the 
pulse to turn off the local oscillator (LOP). Then the switch on of the trans- 
mitter radio frequency is controlled (RFC), which is phase flipped (FLP) between 
0° to 180° for coding. The analogue-digital-converter starts sampling (ADC) 
after RFC off and the TRX, RGT and LOP have opened the receiver and after the 
reset of the Integrator address (RIA) has taken place. At certain range gates a 
calibration signal (CAL) is injected. The whole sequence is repeated after one 
interpulse period when only the phase flip is inverted (for DC-elimination). For 
application of the complementary coding scheme a second double pulse code (con — 
sistlng of the complementary pattern) is transmitted, and this radar cycle is 
repeated N times to yield one radar burst. Additionally, different frequencies 
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(FR.) and pulses for antenna control (ANT) and other purposes can be generated, 
which can change from one to the next Interpulse period or radar burst. In more 
advanced systems multichannel receivers and ADCs are used (e.g. spaced antenna 
parallel processing (e.g. R5TTGER, 1981)). 

The MST radar operation, which Is synchronized and which Is usually started 
and stopped by external clock control (In the host computer), takes place In 
several nested sequences i (1) the radar cycle, l.e. the transmission of one 
code unit of radar pulses with preselected duration and the sampling of the real 
and Imaginary signal at preselected range gates (e.g. 128). (2) the integration 
cycle (burst), l.e. the repetition and coherent Integration of an externally 
selected number of radar cycles. (3) record cycle, l.e. the repetition of a 
preselected number of integration cycles (e.g. 64) to form one total record 
which Is stored in memory. Together with the time and other system parameter 
Information (such as the radar controller program etc.) the data record Is then 
dumped via the computer to file or magnetic tape. 

All the control sequences for the radar system are generated by the radar 
controller, of which a basic block diagram Is shown In Fig. 48. The computer 
loads the Instructions, which pulses should be generated and for how long these 
should be switched on (duration) as well as how many Interpulse periods should 
be coherently Integrated, Into the radar controller memory. The computer can via 
the real-time clock start and stop the radar controller. Further Interlocks, 
either via the computer or other fast hardware devices, Inhibit under certain 
failure conditions radar controller pulses in order to protect the radar system 
from malfunctions or damage if transmitter pulsing would be during the receiving 
phase. 


radar cowtroller/systch synchronizer 



Figure 48. Block diagram of a radar controller or system synchronizer. 









105 


It Is generally found that 8-bit ADCs are sufficiently matched to the line- 
arity range of the receivers (> 60 dB). If strong clutter signals are present 
the capacity must eventually be extended to >12 bit or a range-dependent attenu- 
ation has to be used. Applying an 8-bit ADC and a 1 6— bl t Integrator allows to 
add samples from at least N = 256 radar cycles In one channel. If more samples 
would have to be added, the adder and result memory word length would have to be 
extended. The dump could still be at 16-bit word length if a scaler would be 
used to scale down the result data according to the number of Integrations. The 
preintegration time for 256 integrations would be t* = 2N*Ti PP = 128 ms for Tipp 
= 250 us, corresponding to a maximum unambiguous range of 37.5 km and a maximum 
resolvable radial velocity of 12 ms" 1 (fD»*x = 4 Hz). Applying a 4-bit comple- 
mentary code and a 300-m range resolution would result in a- transmitter duty 
cycle of about 3%. 


MUFFIN FUNCTIONAL BLOCK DIAGRAM 



From Host 
Computer 


Control ler 


Figure 49. Block diagram of a multi-channel finite-impulse response 
filter and integrator (MUFFIN) under development as preprocessor 
for particular coding schemes applied at the EISCAT radars. 

There exist many more sophisticated data preprocessing units at several 
radars and it cannot be in the scope of this lecture to describe those in de- 
tail, as It was not possible to line-out all the different variations of the 
transmitter-receiver-antenna systems. We finally would only like to add the 
block diagram (in Fig. 49) of a multi-channel finite-impulse-response filter and 
integrator which is presently under design at the EISCAT radar facilities and 
shall be used for on-line preprocessing (integration, decoding, correlation 
computation and radar control) of incoherent scatter and coherent scatter data 
(POSTILA, 1988, personal communication). 
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6. AN OVERVIEW ON SOME RADAR SYSTEMS AND THEIR TECHNICAL PARAMETERS 

Descriptions of the early VHF radar systems were given by WOODMAN and 
GUILLEN (1974), GREEN et al. (1975), C2ECH0WSKY et al , (1976), ROTTGER et al. 
(1978) BALSLEY et al. (1980), and design considerations were summarized by 
BALSLEY ( 1978a, b), GAGE and BALSLEY (1978) and BALSLEY and GAGE (1982). More 
recent developments and technical details were described for instance by 
CZECHOWSKY et al . (1984). KATO et al . (1984), ROTTGER (1984), STRAUCH et al . 
(1984), FUKAO et al. (1985a, b), FRISCH et al. (1986), ECKLUND et al. (1988) and 
ROTTGER and LARSEN (1989). 

We present in Table 5 the technical parameters of the MU radar (Middle and 
Upper atmosphere radar) of the Radio Atmospheric Science Center of the Kyoto 
University (from FUKAO et al,, 1985a, b), in Table 6 the basic parameters of a 
405 MHz wind profiler radar (from FRISCH et al., 1986), In Table 7 the primary 
system specifications for an operational wind profiler (from UNISYS, 1987) are 
repeated, which is supposed to be used In meteorological routine applications. 
In Table ' 8 the basic parameters of the EISCAT incoherent scatter radar systems 
operating in the 933 MHz and 224 MHz bands are summarised. We then display in 
Table 9 a most complete list of all the MST, ST and incoherent scatter radars as 
well as some wind profiler systems. Finally we show in Fig. 50 photos of a few 
well known research radar facilities. 


TABLE 5 


Basic Parameters of the HU (Middle and Upper Atmosphere ) Radar 
Operated by the Radio Atnosohere Science Center of the Kyoto University 
(from FUKAO et al., 1985a, bl 


LOCATION 

RADAR SYSTEM 

OPERATIONAL FREQUENCY 
ANTENNA 

Bean width 
Steerability 
Bean directions 
Polarizations 
TRANSMITTER 


Peak power 
Average power 
Bandwidth 

IPP 

RECEIVER 
Bandwidth 
Oynamc range 
IF 

A/D converter 
PULSE COMPRESSION 


Shigaraki. Shiga, Japan 
( 34.85 # N, 138, 10 5 E) 
monostatic pulse radar; 
active phased array system 
46,5 MHz 

circular array of 475 crossed Yagis 
aperture 8330 n 1 1103 i d i ameter ) 

3.6 8 (one way; half power for full array) 

steering is completed in each IPP 

1657; 0 C -30 S off-zenith angle 

linear and circular 

475 solid state amplifiers 

(TR nodules: each with output power 

of 2.4 k* peak and 120 I average) 

1 MM ( maximum) 

50 kM (duty cycle of 51) inaxiiun) 

1.65 MHz (maximum) 

(pulse width: 1-512 us variable! 

400 us to 55 ns (variable) 

1.65 MHz (maximum) 

70 dB 
5 MHz 

12 bits x 8 channels 

binary phase coding up to 32 elements. 

Barker and complementary codes 
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TABLE 8 Tin E1SCAT fitmr ficilitus uan, iMl 
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*) operation discontinued •) in planning or construction phase M ; iesosphere, S = stratosphere, T= troposphere 

nostly ionospheric radar ♦> TX antenna gain in bistatic node (M) = D-region in incoherent scatter node 

» sostly for radar astronoiy $A = spaced antenna capability I = ionosphere-theriosphera incoherent scatter 
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Figure 50. Photos of radar observatories: 
Jlcamarca/Peru, Areci bo/Puerto Rico, 
SOUSY/West Germany, Poker Flat/Alaska, 
HU Radar/Japan, EISCAT /Norway. 
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